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PREFACE 

'' I ""HE Theory of Map Projections has received the attention of 
-*- so many eminent mathematicians and geodesists that it may 
seem somewhat presumptuous to offer anything further on the 
subject. These Papers must therefore prove their own justification. 
They are the result of investigations extending over (with intervals) 
nearly four years, by one who was for twenty years employed on 
the Trigonometrical and General Survey of one of the British Pro- 
tectorates in Tropical Asia, but who had until recently devoted little 
attention to the theory of map projections, being content to follow 
the established methods of the Ordnance Survey and the Survey of 
India. 

The investigations came about in this way. At the beginning 
of iQi6 the writer had offered his services to the Secretary of the 
Royal Geographical Society to assist in the map work being done 
by that Society for the War Office. One of the first tasks given him 
by Mr Hinks was to compute some Tables of rectangular coordinates 
of the special modified Polyconic Projection adopted for the Inter- 
.national Map qn the scale of i/million, for a few sheets to be 
compiled on half that scale. The carrying out of this calculation 
led to an examination of the basis of the formulae on which the 
modification referred to depends. This led to a comparison with 
other Projections and their errors, until, becoming completely 
fascinated with the subject, the writer thinks he can claim to have 
passed under review the principal books and Papers which have 
been published on it, and also to have worked out some results 
which are new so far as he knows, and which he hopes may prove 
useful to practical map-makers and surveyors. These results are 
concerned mostly with the Minimum Error Projections, and although 
the writer cannot claim to have discovered anything of striking 
novelty in either theory or practice, he has been able to extend 
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the method invented by Airy and Clarke more generally, and to 
simplify their results ; so that, out of the almost bewildering number 
of Projections that have been discovered and advocated from time 
to time, those which are practically the best or most useful are 
reduced to comparatively few, the application of which has been 
simplified. One of the most interesting results is to establish in 
their proper place some old and, by later writers, somewhat dis- 
paraged Projections, which seem to have been hit upon by their 
inventors by a sort of geometrical intuition rather than by rigorous 
analysis. 

Owing to the present high cost of printing the Papers are pre- 
sented in as concise a form as possible. They do not profess to be 
a complete treatment of the subject ab initio, but only a statement 
of what are supposed to be new results. Many steps in the mathe- 
matical processes have been left out, and in the last parts it will be 
difficult to follow the argument without reference to the works of 
the other writers mentioned. 

The Tables etc. have been calculated with s-figure logarithms 
and Fuller's slide-rule, and are intended only to be illustrative. In 
the formulae the symbol log means the natural logarithm and the 
angles are supposed to be reckoned in circular measure. Though; 
. the work has been checked where possible by alternative methods, 
there has been no absolutely independent check, and it cannot be 
that some errors have not crept in both in formulae and arithmetical 
work. The writer will therefore welcome corrections from any reader 
who may have the will and opportunity to check his results. 

In the course of these investigations the writer received valuable 
hints as to suitable and fruitful lines of inquiry from Mr A. R. Hinks, 
F.R.S., Secretary, Royal Geographical Society, to whose own little 
book on Map Projections he is indebted for a clear preliminary 
survey of the whole subject. He has also received valuable assistance 
in mathematical difficulties from Professor A. E. H. Love, F.R.S., 
Oxford, Dr T. J. I'A. Bromwich, F.R.S., Cambridge, and Professor 
J. W. Nicholson, F.R.S., King's College, London. He is indebted 
to Captain G. T. McCaw, General Staff, for reading the original 
draft and for making some useful suggestions. 
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CHAPTER I 



THE MINIMUM ERROR ZENITHAL PROJECTIONS 



The Zenithal Projections are a particular case of the Conical Class 
obtained by opening out the cone until it becomes a plane or, which is the 
same thing, by making the constant of the cone equal to unity. They will 
be considered before the general conical class, however, because they are 
the first to which the principle of Minimum Error was applied, and it will 
be simpler to approach the subject from what has already been done. 

Suppose that r is the distance from the centre or pole to any point on 
the map, then if Q is the corresponding angle subtended by this distance 
at the centre of the Earth, supposed spherical, we have as the law of the 
projection r=Rf(&), where R is the radius of the Earth. In the following 
investigations we shall for simplicity assume ^ = unity, but it must be 
remembered that to convert the values of r in any case into actual dis- 
tances, they must be multiplied by the proper value of R on the appropriate 
scale. 

The radial scale of the map is given at once by dr\dQ and the perperi- 
dicular scale by i^/sin 6. 

By choosing different forms of the function r =/{&) we obtain Zenithal 
Projections with different properties of which the following are the principal : 





Name 


r=f(e) 


Radial Scale 
= dr\dd 


Perpendicular Scale 
= r/sin e 


r. 


Equidistant 


e 


I 


cosec 6 


2. 


Orthomorphic 
or Stereographic 


2 tan ^0 


sec^ \6 


sec^ P 


3- 


Equal Area 


2 sin ^d 


cos|^ 


sec ^6 


4- 


Gnomonic 


tan^ 


sec^(9 


sec 6 




Perspective 


ksmO 
h + cos 6* 


k {h cos e+i.) 
(h + cos ey 


k 


5- 


h + cos Q 



It will be observed that the fifth only is multiplied by a constant k, 
commonly called a scale constant, but it will be shewn later that the proper 
evaluation of this constant plays a most important part in the determination 
of the criterion of minimum error, and that a similar constant could be 
applied to the other four with like result, which, so far as the writer is aware, 
has not hitherto been noticed. The second projection in the table has the 
property that the scale in any direction round a point is the same, and the 
third the property that the radial and perpendicular scales at any point are 
reciprocal, so that equality of small areas round any point is preserved. 
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Airy's Zenithal Projection by Balance of Errors. 

Sir George Airy, late Astronomer Royal, conceived the idea of a pro- 
jection which should be a sort of compromise between these two. His 
investigation is published in the Philosophical Magazine for December 
1861. It consists in finding by mathematical analysis the form of the func- 
tion r=f{6) which makes the squares of the two scale errors, integrated 
over the surface of the sphere represented by the map, a minimum. Airy's 
method consists in fact of making the average scale error, or more strictly 
speaking the mean square scale error, of the map a minimum ; this average 
being found by squaring the scale errors in the two directions mutually 
perpendicular, integrating them over that portion of the sphere represented 
by the map, then (implicitly) dividing by the area of the same surface and 
extracting the square root. The result will be the mean square scale error 
of the map, from which the average scale error can be found by multiplying 
by 0-7979.... 

Unless orthomorphism or equality of small areas is wanted for some 
particular purpose, it is difficult to conceive a more rational criterion than 
Airy's for finding the best projection. Another criterion is to make the 
maximum positive scale error equal to the maximum negative scale error, 
in other words to make the maximum scale error as small as possible. 
Since the scale error increases as the square of the distance from the pole 
of projection, which is the centre of the map in the zenithal class, this has 
the effect of diminishing the scale error at the margin of the map at the 
expense of increasing it near the centre and increasing therefore the average. 

Nevertheless, Airy's Projection by Balance of Errors has not been much 
used, and it seems to have attracted more attention on the Continent than 
in this country. The reason, probably, is the apparent complexity of its 
formula, but also, the writer thinks, because its nature has not been fully 
understood or explained in the text books — the more so since Airy himself 
made a slip in his solution which caused the projection to appear under a 
disadvantage in the Tables published with it. Although this slip was dis- 
covered by Clarke and rectified in a paper published in the Philosophical 
Magazine for April 1862, Airy's erroneous tables have unfortunately been 
incorporated in later publications, e.g. A Treatise on Projections by Thomas 
Craig, United States Coast and Geodetic Survey (Washington, 1882), and 
others. 

It also seems to be implied by some writers tliat Airy's is the Minimum 
Error Projection, whereas it is only the Minimum Error Projection of the 
Zenithal Class, which is a particular case of the more general Conical Class. 
The writer has succeeded in applying Airy's method to this more general 
class, and the results derived will be given later. 

Germain gives a full account of Airy's Projection with Clarke's correction 
and in his general review p. 227 says, " La projection de sir Airy, Balance 
of Errors, n'a pas encore ete employee ; sa formule, tres-eompliquee ndces- 
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site des calculs tres-penibles, et les problemes de la sphere ne peuvent 
s'y resoudre facilement." Its qualities and the modification of it called 
Breusing's Projection are also treated by Hammer; but the writer has 
found by far the most complete treatment, with some extensions of Airy's 
and Clarke's method, in Le Projezioni delle Carte Geografiche by Pro- 
fessor Matteo Fiorini (Bologna, 1881), for reference to which he is indebted 
to a pamphlet by J. I. Craig of the Egyptian Survey Department on The 
Theory of Map Projections (Cairo, 1910). 

. Returning to Airy's solution, the two scale errors at any zone at 
angular distance Q from the centre of map are drjdd-i and r/sinfl — i. 
The corresponding element of area on the sphere is 
Rde . 2-a-R sin B = 2-kR'' sin U&, 
and the problem therefore is, omitting the constant multiplier 2-ieB?, to find 
the function r=f{B) which will make the integral 

^=/o lU" ^ (sb" } ''" ^'^^ ^ minimum, 
where /8 is the angular radius at the limit of the map. The corresponding 
area on the sphere is 2TrR^ I sm ^ ^^ = i — cos /8, omitting the constant 

multiplier 2irR^ ; so that the square of the mean square error, which we 
shall call /i^, is given by Ml{i — cos P). By the Calculus of Variations, Airy 
finds 

r^A tan \6 + 2 cot \6 log sec \6 + B cot \d, 

where A and B are constants to be determined by the boundary conditions. 
Since r must = o when 6 = 0, Airy and Clarke put at once B = o. Airy then 
makes the mistake of supposing that the constant A is completely at his 
disposal and he chooses A so that the scale error shall vanish at the centre,. 
i.e. he takes A= i, so that 

r = tan \9 + 2 cot \6 log sec ^6 simply, 
from which his Tables have been calculated. The Calculus of Variations 
exacts however that at each boundary dr/dd=i, and we ought, strictly 
speaking, to suppose the map to extend from a lower limit 6 = a to a higher 
limit ^ = yS, at each of which dr/d6 =1. 

From these two equations we should obtain 



_ /log sec ^^- log sec M . 
V tanU^-tan^Ja J' 



/ tan" ^g log sec |;3 - tan° ^yS log sec |a \ 
V tan^yS-tan^^a )' 

if then our map extends from the centre we put a = o and find 

A=z 2 cot^ 1/8 log sec J/8, B = o, 
which are the values of the constants found by Clarke. If our map were to 
extend over a zone of the Earth's surface from a spherical radius a to a 
spherical radius y8, we should reduce the average scale error by making use 
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of the constant B as calculated above, and the term B cot \&. It is not 
likely, however, that such a map would be of much use in the zenithal 
system, though, as we shall see later, this term is of great importance in the 
zonal maps of the Conical Class. 

Let us now examine the separate terms of Airy's solution. The first, 
A tan \Q, is evidently half the characteristic formula of the Orthomorphic, 
y = 2 tan \Q, or it is the solution of the differential equation drlM = /-/sin Q. 
The second appears more complex, but if we differentiate it we find 

— (2 cot \Q log sec \Q)=\ — cosec^ \6 log sec \B, 

and if we divide it by sin B we find 

2 cot \Q log sec ^5/sin = cosec^ \Q log sec \6, 
so that 2 cot \Q log sec ^^ is half that function which has the property that 
its radial scale error is equal and opposite to its perpendicular scale error; 
It is in fact half the function which solves the differential equation 

{drjM - i) + (r/sin 61 - i) = o, 
that is, r = 4 cot \6 log sec \Q + B cot \B, 

where in the present case B = o. It is therefore to a first approximation 
half a function giving an Equal Area Projection, the characteristic function 
of which is /•=2sinJ5. This is seen more readily from the expansions of 
these two functions, viz., 

4 cot le log sec ie = e- ^^ - ^^e^ -.... 

We see therefore that Airy's Projection is in form practically a mean 
between an Orthomorphic and an Equal Area Projection, but the ortho- 
morphic part must be multiplied by a constant A which has the effect of 
reducing the positive scale error at the outer boundary, at the expense of 
introducing a negative scale error at the centre. The function 

r=4 cot |fl log sec ^6 
does not seem to have been used by itself for a map projection, nor does 
there seem to be any reason for preferring it to the simpler and practically 
identical Equal Area function /- = 2 sin ^6. It lends itself more readily, 
however, to the mathematical operations necessary to determine the con- 
stants of the Minimum Error Projection and is a useful guide when dealing 
with the absolute Equal Area function in both the Zenithal and Conical 
Classes. 

By substituting the value of r as found above in the expression for M 
and integrating, we can find the value of M and hence of /a in terms of /3. 
We shall find 

^-.2sinH;8-^('°gffj^^\ 
^'^ tan^/S 

2 71^1/ a\ 4 (log sec i;8)^ sec is 
ft' = M (i-cosB) = i-3_!^-2 — . -,' si:. 
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Clarke in his Paper quoted above does not appear to have taken this stepi 
but he evaluates the value of M for the particular case of ;8= 113° 30' by 
mechanical integration, getting the value J/= 0-1569. The above formula 
gives o'i5688, which checks Clarke's value exactly. 

Airy's Solution expanded. 

It is instructive and more useful for the purpose of comparison with 
other solutions to expand Airy's solution as follows : 

We have r = A tan -6 + 2 cot - B log sec - Q. 

2 2 ° 2 

Now tani|9 = ie + — 63+_J_5)5+_i7_g7^ 

2 2 24 240 40320 

2 cot - (9 log sec - 61 = - e - 4 ^' - -^ 6' ^— &•-..., 

2 2 2 48 960 17920 

sin 61= 0-\ (PA-^^ff' ^ 6'^. _ 

6 120 5040 

We know that the first term of the constant A expanded in terms of )8 is i, 

and we shall save some trouble by writing A='l +b, where b only is to be 

found. We can then add the first two expansions together and write 

48 320 101,280 2 V 12 120 / 

thus 

drldd=i + ^e'+^ 6"+-^^ e'+...+-b (1+ 1$'+ — 6"...), 
16 64 23040 2 \ 4 24 / 

Wsin^=i+4^'+ — ^' + — ^ ^' +-^ (1 + - 6^+— 0' ); 

' 16 192 23040 2 V 4 24 / 

therefore 

M= r {{drjdO -■Lf + (r/sin 6 -if] sin 6 dd - 
Jo 

Jo (.256 \ 4 1440 / 

16 V 4 1440 /\ 4 24 /4\ 4 24/ 

256 V'' 12 1440 / 16 V'' 12 1440 /\ 4 24 / 

jo 1256 V 3 9 / 16 V 60 y 4 v 3 120/J 

= I^V3 "6 45/ 16^ 6 ^ 480^^ 4^6 12oJ 
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Putting the differential coeflficient of this with respect to b equal to zero we 
find for the value of b to make this expression a minimum 



8 11520' 

and consequently yi = i — '^ + — - — , 

^ ' 8 11520 

which will be found to agree with the expansion of 2 cot^ \^ log sec \^ as 
found above. 

Substituting this value of b in the expression for M we find 

M = -^ + -^^— -^ --^^ — 



768 1536 11520 

256 1536 11520 ■■■ 

/8° 
= -^ to ■? terms expanded, 

384 

^ I ^ \ 
dividing M hy i — C0S)8 = — (i — ^— + ...) we find ft? and consequently 

02 / B' \ B' 
a = — f- — I + ^— ... I = ^— nearly. If we substitute the value found for b 

^7192 V 24 y 14 •' 

in the original expansions we find for the scale values at the centre, 

I B^ 

where 6 = 0, {drjdO^ = (r/sin 6),,= i + -b = i-^ + ...; and for the limiting 

radius and scale values thereat 

^"'^V 24 480 8960 "V' 
{driae)p = I, 

B^ B* 
(./sine), = x + |+|+.... 

Supposing instead of using the value of the constant A which makes 
the mean square error a minimum, we had used ^i = i, as did Airy, then 
we shall find the values of the corresponding quantities by making b=o 
in our formulae ; thus (drldB\ = (r/sin ^)o = i, 

48 320 ■ 
/3^ /8* 



-K-S^l^-)^ 



^'^i''^^- ^^z6 64- 

(/-/sin (9), = I + 5^% M! , 
16 192 



so that both the maximum and mean scale errors are about 50 per cent, 
higher. Airy's solution with A = i does not even give the minimum error 
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solution with the stipulation that the scales shall be true at the centre, as will 
be proved later, which shews the importance of correctly evaluating the 
constant A. '' 

To find for what value of 6 in the Minimum Error solution the perpen- 

dicular scale error vanishes put Wsin 6 = and we find ~ . . . — ^ = o or 

16 16 

B = —j- /3 nearly. 

It might appear that the expansion for M contains no terms beyond /S"; 
there is, however, a term in /3^^ as follows : 

384 V 3024° / 

Breusing's Projection. 

As we have already remarked, Breusing conceived a projection which 
is the geometrical mean between the Orthomorphic and Equal Area, viz., 
r= 2 y^^sin ^6 tan ^0), and Hammer has pointed out how nearly this agrees 
with Airy's solution taking A = 1. The writer finds that whether the 
arithmetic, geometric or harmonic mean of these two functions be taken, 
we obtain a projection agreeing with Airy's to 2 terms in the expansion of 
M, and to at least that order in the other quantities. The harmonic mean 
leads to a much simpler function than either of the others and we shall 
investigate it accordingly. 

The Breusing (Harmonic) Minimurn Error Projection. 

The harmonic mean of 2 tan ^6 and 2 sin ^6 is 4 tan ^6. Consequently 

write r= \A tan \Q and determine the best value of A as before. We shall 

find 

4 - si"'i<^ M- A .;in= A/? - i^sin^i^ 

logsec^;8-tan^i;8' 4 sm ^P log sec ^/S - tan4^ ' 

To shew how nearly this projection agrees with Airy's we shall find in 
applying it to the extreme case taken by Clarke, viz., /3=ii3°3o', the 
limiting radius comes out as i'578i as against f576o, and the value of M 
as 0-1579 as against 0-1569, i.e. its average scale error is only i part in 
3000 greater. 

By expanding it in the same way as Airy's we shall find ' 

256 



3 


12 5760 


— I — 


^. 


/8^ 




12 


120 



) 384V 3840^ 



For the scale values at the centre we get 



{drldO)o = {rjsm 0\ = A = i-^ 7g , differing in the 3rd term only. 
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At the limiting radius we find 

''^ = ^(' - £■- J^- 8^^-)' ^S""'"'"^ '° 4 terms with Airy's, 

(r/sine)^= 1+^+^... „ „ 3 

This much simpler formula thus leads to almost identically the same 
result as Airy's, and when we remember that in the extreme case of 
/3= 1 13° 30', a map covering more than two-thirds the surface of the globe, 
the maximum error is actually 72 °/„, it certainly does not seem worth while 
using the more complex formula of Airy to reduce the mean error by only 
i/30ooth part. 



The Orthomorphic Minimum Error Projection. 

In this we have r=2A tan \Q and we shall find that A has exactly the 
same value as in Airy's and M double that value, viz., 

A = 2 cot^ JyS log sec ^/3 = I - ^' + ^' 



11520 



tan^^^ 192' '^ 796 V 24/ 

The scale value at the centre is 

(dr/dd), = (r/sin e\=.A=j-^ + -^ . .., 

8 11520 

and the limiting radius and scale value thereat have the same values as in 
Airy's. 

The Equal Area Minimum Error Projection. 

In this we have r=2A sin ^6 and on integrating we shall find 

^ _ |(i - cos^ i^)+i- cos i/3 p' 

i(i-cos4/8) + logseci^~^ I28'"' 

Hi-cos*^/?) + logsecJ;8 192 V 16 768 T^) 

192V 16 192/' 
the term -— ^ bemg the diminution in M caused by the use of the con- 
stant A. This system is therefore capable of practically no improvement 
by the use of a scale constant. 
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The scale values at the centre are ^4 = i - ^-— .... and at the outer limit 

128 

of the map we have 



'■? = 


..(■ 


24 


7/8^ 
960/' 


{drldO),- 


I 


13- 
8 


192 


rjsm 6)^ = 


= I 


8 192 



The use of a scale constant with 2 sin ^6 is in fact not appropriate to 
this Projection, for, as we see, it destroys the characteristic property of 
reciprocity of scale values, at any rate after the first 2 terms. The same 
xemark applies to its prototype r = 4 cot ^6 log sec |(9 for which also 

-M= . It will be shewn later that for this we should write the full 

102 

formula r = 4 cot |6 log sec ^d + 2B cot \A, and after determining B so as to 
make M a minimum for a zone from 6 = a to 6 = /8, we put a = o, which 
amounts to making B = o originally. 



The Equidistant Minimum Error Projection. 

Clarke in his Paper on Airy's Projection published in the Philosophical 
Magazine for April 1862 compares it with two others. One, the Perspective 
Projection of James, has received considerable notice but the other, the 
Equidistant, or Equal Radial as Clarke calls it, is seldom mentioned, 
although as he shews when drawn with the same limiting radius as Airy's 
its M is to Airy's as 1047 : loo'o. Clarke does not explain very clearly 
-why he draws the limiting radius equal to Airy's, but what he does has the 
effect practically of writing the equation of the Equidistant as r = AQ and 
then of determining A so that the mean square scale error shall be a 
minimum. 

On applying our method we find 

i-cos;8 + |^^ _ /6^ 7/8^ n^' 

'^/,, /ij/i 24 2160 362,880'"' 



Jo 



I-C0S/3+ Vcosec^^^ ^^28^ 75 "8,125 ) 
.'0 
For p= 113° 30' this gives iI/=o'i634 as against 0-1643 found by Clarke 
by mechanical integration, for the case r^ the same as in Airy's solution. 
The scale values at the centre of the map are 

A = r-^-.^^' 



24 2160 
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and at the outer boundary we find 

Va = fi ('i - — - -^V- • ■• V agreeing with Airy's to 2 terms, 
\ 24 2160 / 

{drlde)^= 1-^ ^... „ „ I term, 

' "^ 24 2160 

■ EP B* 
' ' 8 108 

The value of the mean square error /j. 

12 

This exceeds the value in Airy's in the ratio Jig2 : Ji'j2-8 = 105 : 100- 

In the writer's opinion this slight increase in the value of the average scale 

error is far out-weighed by the valuable property of uniform radial scale- 

and the simplicity of the formula of the Equidistant Zenithal Projection. 

fl* kB^ fi" 

If we make ^ = i we shall find M= -~ as against ^7^ = 7 in Airy's 

216 " 768 iS3'6 -^ 

formula. Therefore Airy's Projection as exhibited by the author and other 

writers is not per se so good as the Equidistant in the ratio 

^216 : Jj53-6= 119 ; 100. 

The merits of this projection in comparison with others are described 
in a Paper by General de Coatpont in the Bulletin de la Socikt'e de Geo- 
graphie for 1894, but without the use of the scale constant. Its invention, 
is attributed to Guillaume Postel in 1581. 



Modification of Minimum Error Projections by making 
total area of map true. 

It will be observed that in all the Minimum Error Projections considered 

the formulae for the limiting radius r^ when expanded agree to the first 

• / fi^ \ . . . 

two terms, viz., /■p = /S(i — ^— ...1. Now in the Equal Area Projection 

ra = B { I — ^— H — . . . I . It therefore appears that we shall sacrifice 

'^ "^ \ 24 1920 J 

very little in the value of the average scale error if we choose our scale 

constant so that equality of total area is absolutely ensured, and this is- 

a valuable quality in any map. 

In some cases, as we shall shew, the formula for the quantity A is 
simplified. In every case considered the expansion for M is altered only 
in the term /3^° and upwards. 
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(i) Airy's Projection with total area true. 
Here " rp = A t&n^ 13+ 2 cot ^/S log sec |yS ; 

therefore to make total area true equate this to 2 sin |^/3, whence 

A = 2cosAfi-2cotUSlogsecAfi= i-^ + — ^ .... 

Substituting this we find 



384 \ 384/ ii92\ 24 

/■p = ;8 f I - -^ + -^V 2 sin 1/3, 
'^ ' V 24 1920/ "^ 

(drlde)p = i + ^; Wsine)^= I + 1+ |g+ ... =seciA 

(ii) Breusing (Harmonic) Projection with total area true. 
Here r^ = 4A tan Jj8 = 2 sin |/3 = 4 sin l^ cos ^/8, 

whence ^ =^°^^i^ = ^ "7^ + ^ - ^il6^+ -' 



24' 



(<fr/<f^)3 = I ; (r/sin e)p = I + ^ + |g . . . = sec J/3. 

(iii) Orthomorphic Projection with total area true. 
Here rp = 2 A tan \P=2 sin \ (3, 

therefore • ^ = cos J^= i - g- + p^- ^g^^+ ..., 

whence ^ ^=^ (^ + ^s) ' '^ = TP V ^ M -J ' 

(^r/rf^)p = (r/sin 6)^ = I + ^ + ^ ... = sec 1/3. 

(iv) The Equal Area Projection of course remains unchanged and 

we have 

rp = 2sin^p, A = i, 

(<fr/./^)p=i-|+ J^...=cosJA 
(/-/sine)3=i+f +|g... = seci^. 



12 
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(v) Equidistant Projection with total area true. 
Here r^ = /4/3 = 2 sin Jj8, 



whence ^ = ^sin|;8=i 



24 1920 



1728 V 75 10,080,000/' '^ V V864/'^ V 200 



7 



(drjde)p = i-tL 



^ P' 



24 1920 



■■=-gSin|/3, 



(r/sine)^=i + |'+|g +...=sec|/3. 
The perpendicular scale error vanishes when = ^/3 nearly. 



The Minimum Error Perspective Projection. 

The geometrical meaning of the Perspective Projection is as follows : 
Let the circle in Fig. i with centre O represent a section through the 

P 




Fig. I. 



centre of the Earth, and draw a line through O and a point M on the 
surface which corresponds with the centre of the map. Then choose a 
point F on the diameter MO, prolonged if necessary, from which points 
on the surface of the Earth can be projected by straight lines on to a plane 
of projection pp which is at right angles to the line MOP. Let A be such 
a point and A' its projection, then the distance JVA' is r in our formula, 
and the angle JVOA is 6 ; denote the distance OF by h and the distance 
NF by k. The radius of the Earth OA is taken as unity. Draw AQ 
perpendicular to MOP and let N be the point where the line MOP cuts 
the plane//. Then 

4IL^^ .„. ^'iv-.^.m-y^AQ ^sin^ 

AQ 



NF 

j-p and A'N=r-- 



QF OF+QO h + cosO' 
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which is the characteristic formula of this class. It should be noted that 
both Germain and Thomas Craig make the mistake of saying that k is the 
distance of the plane of projection from the centre, or rather Craig says, 
"...making h the distance of the eye from the centre and k the distance of 
the plane of projection from the same point...," whereas k is the distance 
of the plane of projection from the eye (i.e. from point F). 

If the point F is taken at F' on the surface of the Earth and the plane 
// tangent at M on the opposite surface, then h=i and k = 2, so that 

2 sin ^ , . 

r — J, = 2 tan A d, 

I + cos 6( ^ 

and we get the ordinary Orthomorphic or Stereographic Projection. 

If we take h = o and ^= i we get r= tan 6 the Gnomonic Projection. 

If we take ^ = i + A so that the plane of projection becomes tangent at 
M, and then make h increase without limit so that the projecting rays 
become parallel, we get ^= sin Q, the Orthographic Projection. 

Several Perspective Projections have been put forward from time to 
time with the plane of projection tangent at M, or taking k=\^h, with 
values of h ranging from i to 273, e.g. those of Parent, Lowry, La Hire, etc. 

General Sir H. James, Director General of the Ordnance Survey, in 
seeking for the best values of h and k for a map covering two-thirds of the 
sphere, or y8=ii3°3o', either by intuition or luck hit upon the value 
// = I '5 and he took the plane pp at the same distance from the centre of 
the Earth as the plane of the tropics, that is, he made 

jVO = sin 23° 3o' = o'399..., 

and therefore k = NO+ 0F= 1-5 + 0399 = 1-899.... 

Fiorini shews (pp. 213-4) that with these values of h and k James' Projec- 
tion gives 

rp= 1-5812, M=o-i']20, (drjdd)^ = o-62g2, (r/sin ^)^ = 1-7243. 

Clarke, by taking A so that r^ = 1-5760, the same as in Airy's Projection, 
finds ^=0-1718, which differs very little from the value with James' values 
of k and A; but on finding by integration the values of k and i which 
make M a minimum Clarke obtains A= 1-36763, A= 1-66261, r^ = 1-5737, 
j1/= 0-1626, {drldO)^ = o-8oS3, (;-/sin 6)^ = 1-7160. 

Clarke on integrating the value of M brings it into the form 



4 


sin^/3- 


A' 






,2_ 


<i 


-^y 


3 


il! 



where 

iP={(..,)logf-(.--.)(i-i)}, 
and where m = h+i, n = A + cosp. 
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In the article on Mathematical Geography in the Encyclopaedia Britannica 
(Ninth Edition) Clarke puts these expressions into a more concise form, 
but the writer has found the above expressions easier to compute with a 
table of Squares, Cubes and Reciprocals like Barlow's. The differential 
coefficient of the above expression with respect to h is unmanageable, and 
Clarke then proceeds to find by interpolation the value of h which makes 
E^IA a maximum. The writer in checking the calculated values of log B^IA 
in Clarke's example discovered a slight error in the value for h = \ -39 
which should be 0-4207 12 instead of 0-420665. This correction leads to 
uniformity in the second differences and gives h - i -368 for the value of h 
corresponding with the maximum. 

Perspective Projection expanded. 

By expanding the formula r = -^ from the outset the writer has 

•^ ^ h + cos Q 

succeeded in finding expressions from which to compute the best values of 

h and k directly, as follows : 

Write k = {h+\)(\-V), so that b represents the distance from plane 



tangent at point M to actual plane of projection pp. Also write c = 

and expand cos Q. Then we shall find 

whence 

2 V 3/ 4 V 2 30/ 8 \ 3 10 630/ 

... e^ eW i\ 6« / , r I \ 

Wsm d=\-v—c+ — c\c--\+-=rC\^ + — )... 

' 2 4 \ 6/ 8 \ 3 90/ 

Forming the quantity 

M= I {{dr/de -iy + (r/sin 6 -if} sin edO 
Jo 

in the usual way and integrating, we finally arrive at 

+ M.Y3S^_S9£'^S4S^_9f^_7_\ 1 
S V 2 3 72 8 144/"/ 
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-'f{^('-j)^f(— f) 

(. 2 V 3/ 3 \ 4 4 120/ J 

Differentiating this with respect to b we find the value of b to make M 
a minimum is 

Substituting this we find 

+ «4 /^ _ S9£' ^ r°9^' _ ^ + _7_\ I 

V 2 5 24 40 240/ ■■■/ 

V 6 24 40 480 20/ '"J 

+ 54/1' 3£!_?3^ 83£_J_\ 1 

"where the first two Hnes represent the value of M when b = o, or the pro- 
jection is made on the tangent plane as supposed by Tissot and many 
writers ; the second two lines represent the decrease in M due to moving 
the plane to its best position b units towards the centre of Earth ; and 
the third two lines represent this best value of M. We can now find the 
best value of c on either supposition, by differentiating the first or third lines 
with respect to c, and equating the result to zero. 

Taking the general case first we get for a minimum 

solving this by successive approximations we find 

8 256 61440 8 \ 32 23040 / 

Substituting this value of c we find 

384 V 3°72 / 
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which differs from the same quantity in Airy's Projection by the term in 
j8i« only. We shall also find 

= I -^ - -5i-- '-—- .... agreeing with Airy's to 2 terms, 

16 1536 24576 

^^ '^V 24 480 8960' 

{dr/d0)^=i-^... „ „ 2 „ 

IP B* 
(r/sine)p= I + ^ +^ ... „ „ 3 I. 

and we also find 

c 3 V 20 28800 / 

We shall now check these formulae by computing the values of /t, k 
and M for the cases considered by Clarke and Fiorini, viz., /3= 113° 30'' 
and;8=9o°. We find 

^=ii3°3o' /3=9o° 

h= I '367 as against i'368 /%= 1-471 as against 1-470 

k=\-(i(io „ i'663 ^ = 2-040 „ 2'o38 

M=o'\(ii7, „ o'i626 J/=o"o395i ',, o'°39S9 

The above expansions for h and k are therefore sufficiently exact to- 
enable us to calculate the constants of the Minimum Error Perspective 
Projection for maps of the largest extent, and to calculate them directly 
instead of by the longer and more troublesome method hitherto employed. 

Perspective Projection on tangent plane. 

If the projection is made on the tangent plane we have 

P i/,~^ <:, I ,\ , m(,^^ 43*^ 



\ 2 5 24 40 24o/""J' 

whence we shall find for the minimum 

10 \ 200 120,000 / 240 \ 200 / 

3 \ 200 / 3 \ 20 / 

k = h + \=— U -|^...y Also {drlde\ = {rls.m6\ = i; 
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These properties differ considerably from those of Clarke's Minimum Error 
Perspective Projection. The value of M is however smaller than for any 
other Zenithal Projection considered, with scale values true at centre. With 

yS = - we find h= 2-04, which does not differ much from the value 2 '105 of 

the Third Projection of Parent. 

From the general formulae just investigated we can find the properties 
of the other well-known Perspective Projections on both the tangent and 
Minimum Error planes. For instance 

h = o, or c = I gives the Gnomonic, ;• = ( i - 15) tan Q, 
h='X),oxc = o „ Orthographic, r= {i—b)sin0, 
h=i, or c = ^ „ Stereographic, ^ = 2 (i — 1^) tan ^6, 
the appropriate value of b being found by substituting the proper value of 
c in the above expression for b, if the projection is desired on the minimum 
error plane. It may be noticed that if we take c=\, we obtain a Perspec- 
tive Projection which is approximately Equidistant, for it agrees with the 
latter to terms in 6^ and /3^ inclusive. And as the areal scale is to two 

terms expanded dr/d0.r/smO= i -1 — (4^-1)... we can get a Perspective 

Projection which is nearly Equal Area by taking c= \, or A = ^. In this 
case i5 = o to a first approximation, a result pointed out by J. I. Craig, p. 58. 
Though it has all these inherent properties the Perspective Projection 
has, in the writer's opinion, been over-rated in the text-books written for 
practical cartographers. General Sir H. James extols it above Airy's because 
it is a real, or purely geometrical, projection. But as we have seen, it is in 
its general form awkward to compute both for plotting and for finding the 
scale errors, and it differs inappreciably from Minimum Error Projections 
of simpler form such as the Breusing (Harmonic) or the Equidistant. 

The Minimum Error Zenithal Projection with scale 
true at centre. 

The fact that Airy's solution with the constant A=i does not give the 
Minimum Error Zenithal Projection, subject to the condition that the scale 
shall be true at the centre, caused the writer to try and find why the 
ordinary process of the Calculus of Variations failed in this instance. After 
solving the Perspective Projection with this stipulation it seemed clear that 

the requisite function when expanded would be of the form r=6 — ~6'^..., 

because the Perspective Projection agreed so closely with Airy's in the 
general case. Also if we assume r= 0+ a8^..., where a is a constant to be 

determined by making Ma. minimum, we find J/=^ f loa^ + - + — j , 
whence a = - 7- and M= , agreeing so far with the Perspective solu- 



1 8 MINIMUM ERROR 

tion. On substituting this solution in the fundamental formula of the 
Calculus of Variations 

8 U= {hM in our case) = \ KhrdO + ff^-If„ 
Jo 

dP 
where K=N-^+...; II^ = SrpiFf,- ...) ; H, = K{P,...), 

the equation 8 17= o was found to hold, though the usual separate conditions 
K=o and H^—If„ = o did not hold. 

The writer was fortunate enough to interest Professor J. W. Nicholson, 
F.R.S., iri this question, and with his assistance the fundamental equation 
of the Calculus of Variations was, at any rate partially, reconciled with the 
above indirect solution in the following way. We have 

M= ( {{dr/dO- if + (r/sin d- if] sin Odd, 
Jo 

and in the usual notation of that Calculus 

SM=S [ VdO = ( (N?,r + FSp) dO = o. 
Jo 70 

Now assume r can be expanded in the form r = a6 + d6^ + c6^..., where 

a, b, c ... are constants to be determined so as to fulfil the conditions 

prescribed. Then since drldQ = i when ^ = o we have a = i and 

Now 

iV=aF/3/-=2(^-/sin6i-i), and P=dVldp= 2 sind (dr/de - i), 
also 

p = drldd = I + 3^6^+5^61^.., r/sin 6 = 1 + (b + ^ O'' + (c + ^-^+ -^-^ 6\ 
and we can write ^r=U6'^ + M\.., 8/ = 38^6^+ 58^6*... . Substituting these 

= .[8.{(ro..i)f.(r.-|.^)f}.8.(r..l)f...]. . 
For this to vanish for every variation of b we must have 

6 4 \ 3 360/ 

and if we can say that b must be independent of c, then c = o and 



60 320 60 \ 16/ ' 
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and we shall find 

o \4o 1600/ 240 \ 40 / 
This is certainly smaller than the value of M found in the Perspective 

Projection with scales true at centre, viz., -^—- ( i + ^^ ... 1, but the 

240 \ 200 / 

evaluation of c and higher coefficients, seems to leave something to be 

desired as a matter of principle. So far as the cartographer is concerned 

the formula r=d — -^6^ will be accurate enough for all practical purposes, 

or if desired the next term in the formula of the Perspective Minimum 
Error Projection with scale true at centre may be used, viz., 

^ 6^ 2i6'/32 ffi 

r = d--^ + — .... 

60 4000 150 

This must give a result very close to the absolute minimum of average 
error. The theoretical difficulty is mentioned in case it may attract the 
notice of mathematicians interested in the Calculus of Variations, which 
would appear to require development in those cases where conditions are 
imposed on the derivatives of the function sought, at the limits. 

Method of evaluating the average error. 

There is another theoretical point which deserves notice. It will have 

been remarked that Airy and Clarke integrate the squares of the scale 

errors over the surface of the Earth corresponding to the map concerned. 

It appeared to the writer to be more natural to find the average error by 

integrating this quantity over the surface of the map and then to divide by 

the area of the map before evaluating the minimum. It appears that when 

once the form of the projection has been fixed it is as easy to use one 

method as the other. For example taking the Equidistant Projection 

r = Ad = {i+b)6, the scale errors will be expressed as before but the 

dr 
element of zone on the map will be expressed by 2irrdr=2Trr-j2.dd, so 

du 

that we multiply the squares of scale errors by rdrjdO instead of sin 6. In 

the present case rdrjdO = A'^O. Our integral M is however divided by 

[^ dr R^ 

277 r—.de = 2irA'' ^ , 
j^ du 2 

so that the A^ in numerator and denominator cancel out and we simply 
have to multiply by instead of sin before integrating and differentiating 
for b. By using the expansion method we obtain in fact 

L216 \ 40 14000/ 12 \ 90 1680/ 

'^ \ 12 90/J' 
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and i/=^^rJ;(..MM^)--^;(rH-Sf! ,.il|^^)l 
L216 \ 40 14000/ 576 \ 180 226,800/J 

Li728\ 40 14000/ 1728 \ 180 226,800/J 

-^•i728(,'+ 36o-> 
and since the area of the map = we have 

''=y(A)''^(-'£-) 

found by integrating over the Earth. Integration over the map therefore 

gives a result for ft, greater in the proportion of i + /S760 : i in this case. 

^ 7fl^ 
The maximum scale at the boundary is however i + -x- + ^7- as agamst 

I + ^ + -^ in the other method, i.e. slightly smaller. Application of this 

method to the other fixed functions leads to similar changes in the third 
term of the scale values, but when we come to apply it to the problem of 
finding the minimum error function the question becomes much more 
difficult, as the minimum value of the quotient of two integrals is required. 
With some valuable assistance from Prof. A. E. H. Love, F.R.S., the writer 
was able to put this problem into the form of a differential equation, much 
more complex than Airy's, and an equation involving values at the limits. 
It was not found possible to solve the differential equation in terms of 
known functions, but it was solved in series for 3 terms and gave a result 
agreeing with Airy's solution to 2 terms. The question of finding the 
average error by integrating over the map is therefore not of practical 
importance to the cartographer, but it is mentioned here as a point of 
theoretical interest, and also to reassure those who may be inclined to 
object to solutions differing from the absolute minimum error solutions in 
the third term of the expansion of the functions, by the possibility that the 
validity of those third terms may itself be a matter of argument. 

Summary of results. 

The Minimum Error Zenithal Projections have been thoroughly examined 
by deducing and comparing the formulae for their mean square, or average, 
scale errors. It is shewn that for maps of the largest extent the Airy, 
Breusing and Clarke Perspective Projections are almost identical. It is 
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shewn that the minimum error condition differs only in terms of a negligible 
order from the projection obtained by making the total area of the map true. 

It is also shewn that the Equidistant Projection with total area true has 
an average error only about 5 °/^ greater than the average error of Airy's, 
while the maximum error is the same in both. This small defect is far 
outbalanced by the simplicity of its formula and this Projection is therefore 
recommended to cartographers as the best zenithal projection in all cases, 
except where orthomorphism or equality of area or some other special 
property is required. 

The cases where it is desirable to use Zenithal Projections rather than 
Conical (of which Zenithal form a special case) will be considered later. 



CHAPTER II 

THE MINIMUM ERROR CONICAL PROJECTIONS 

We have already stated that the Zenithal Projections are only a par- 
ticular case of the Conical Class, in which the map instead of being developed 
on a plane is developed on the surface of a cone having as axis the diameter 
of the sphere which passes through the centre of the zenithal map. This 
cone is then supposed to be laid out flat. The only difference is that the 
projection of a small circle on the sphere, perpendicular to the axis, instead 
of being a complete circle of circumference 2-Kr is a portion of a complete 
circle 2Ti:nr, where n is the sine of half the angle at the vertex of the cone 
and is called the constant of the cone. The radial scale is as before drjdO, 
but the perpendicular scale is now nrjsva 0. 

The expression for M becomes 

M= r Fde = r {{dr/d6 -if + (nr/sin O-iY} sin Odd, 

Ja Ja 

and the expression for the corresponding area on the sphere is as before 

(cos a — cos p), leaving the constant multiplier 2TrJ?^ in each. 

In the usual notation of the Calculus of Variations 

JV= d V/Sr= 2n {nr/sin 6 - i), I'=dV/dp = 2 sin 6 (dr/dd - i), 

dF/dd = 2 cos e {dr/d0 - i) + 2 sin 0d'r/de\ 

so that the differential equation iV- dPjdd = o becomes 

d/T df 

sin^ 6 ^Ts + sin 6 cos 6-n^- n^r = sin 6 cos 6 — « sin 6. 
d(r da 

The complete solution of this equation is 

r = A (tan I6f + (cot J 6)" [(tan ^Ofde + B (cot 1^)^ 

which at once reduces to Airy's solution for the Zenithal Class when we 

make n = \, for I tan J 6^5 = 2 log sec \B. The constants A and B are 

determined from the conditions that P=o both when Q = a. and ^ = /3, 
therefore dr/dQ = i when & = a and ^ = yS, or the radial scale is true at both 
boundaries of the map. These two equations lead to the following values 
of the constants : 

[ (tan|(9)»(j 



A = 



(tan i/8)^" - (tan |a)« ' 

(3 



B = 



(tanla)^" f (tan J(9)" (^^ - (tan J/8)2" ('(tan |e)V6l 
h Jo 



(tan |/i)» - (tan ^a)''"' 
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and the limiting radii of the map will be 

2 (tan 1/3)" r (tan \BYd6 2 (tan Ja)" T (tan \efdB 

^^ " ~(t"airpp« - (tan ia)^'' ' ''"" (tan ^/3f " - (tan ^a)» ' 
so that the depth of the map r^ — /■„ 

2 { (ta.n^Ofd0 



! r (tan ^0)''d9 . [ 



(tani/3r-(tan^ar 
(tan 1/3)2" - (tan ja)^" 



(tan^;8)''+(tania)»- 



The only trouble so far is that the integral I (tan \6Yd9 cannot be 

expressed in known functions except when n= 1. It can be evaluated for 
any given value of n in series or by mechanical integration. The calcu- 
lations may appear rather formidable at first sight but they would not in 
any given case occupy an expert computer more than a day or two, which 
is a small matter compared with the total time occupied in the compilation 
of a map. We shall shew, however, that by means of expansions and 
simplifications even this labour can be avoided, at the expense of almost 
imperceptible increments of the average scale error. 

We have still the constant n at our disposal and in order to find the 
value of this which makes M a minimum, it will be convenient to consider 
the Orthomorphic and Equal Area parts of our solution separately. 



Minimum Error Conical Orthomorphic Projection. 

The general formula for this is r = 2 A (tan 1^)", where 

A=i^ (tan \6fdei{{tan IPf- (tan Ja)»}. 

For simplicity we shall write henceforth 

/ for I {ta.n^efde. 

Substituting this in the expression for M formed in the usual way, we shall 
find after some transformations 

M= 2 Fcos a-cos p- 2n( ry /{(tan ^yS)^"'- (tan ^a)=''}1 . 

Writing \{P + o)=X ^"^ h{P~°) = ^> we have 

cos a — cos yS = 2 sin X sin 8 = 2 sin X • S ( I — 7- + . . . j . 

[^ . 
We can also expand the integral / m terms of the middle value of the 

function and the limits by the following general formula 

Ijid) d6 = 28 {/(x) + 1/" (x) + —/- (x) • • •} • 
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In our case 

/(X) = (tan ixr ; /" (x) = ""-^^^i^ (tan W : 

Mil X 

yiv / \ = «' - 6i^' COS X + 7^' cos^ X + 4^' - 5^ cos X - i^ CQS° X ,^^^ :^^n_ 

sin*x 
Thus 

>»(/;)-=s.a.(,„w..{..fJ"^5^ 

S^ i6i7^ — 56/z^ cos X + 52!^° cos^ X + ^4^'^ ~ 3°^^ ^os x — 6^ cos'' x 
360 ' sin* X 

We also have 

(tan ^(Sr - (tan ^aT = 2n I (tan ^(9)* cosec 6^^, 

and expanding this in the same manner we get 
2n (tan Jx)'" f ^^ 4^^^ "^ ^'^ ^os x + i - cos^ x 



sin X I. & sin^ x 



^lxr/i+? 

inx 1 6" 

i6«''-8o«'cosx+ ioo«^cos^x+4°'^'^~7°''cosx— 3o«cos^X -> 



8* + i8cos^X+cos*x+5 



120 sm'x 

Dividing out we find 

2« (£)y{(tan Ipr - (tan Ja)^"} 

„ . ( h'^ 2n^ - An cos X + I + cos^ x 

= 28 sm X i I - 7- • r^^ 

y 6 sm^X 

48;^''— 96;z^X + 52'«^cos^X~3^'^^+8°'^cosx-i6«cos''x -, 

, 8* -34cos^x + 7cos*x-5 l 

360 ■ sin*x "J ■ 

Diiferentiating the last with respect to n to find the turning value we 
get the equation 

8^ 4« — 4 cos X 
6 sin^ X 

8* iga;?'*— 288«^cosx+io4«cos^X~72'' + 8ocosx— i6cos^X 
360' sin*x ~ ' 

whence n = cos x to a first approximation, and inserting this value in the 
second term we find 

i cos X - 8 cos' X _ f 8^ 



6' » cos X - » cos" X ( 6' \ 

« = cos X + • ^Vrj = cos X I + — ■ ■ • I • 

240 sm^x ^ 30 / 

When we insert this value of n in the expression for Af we find 



M-- 



8' 8* 8^ 8* 1 48^ sin X 



:48sinxii-^+ — -...- I +- + — ...^ = 
[ 6 120 6 72 J 



45 
and dividing by the area, fi^ = -^8*. 
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When the map extends to the centre or when a = o we proceed as 
follows. We have 

M=2 \i-co&P-2n(\ \ (cotlpyA . 

Now I = 2 i + . . . L 

;» \i+n i + n 5+« J' 

where / = tan J/8 ; therefore 
« (l^^y (cot i;8r 

'^'^\(T.+nf (i + «)(3 + «)"^ \(i + «)(5+«) "''(3 + «)V '"/ 

= 4^Ho^T^~(i + «)(3 + «)^^((i+«)(s + «) ^ (JT^V ■■■} ■ 

Differentiating the part within the brackets with respect to n to find the 
turning value we get 

from which by successive approximation we find 
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25^^ _. ^%4^\..= j_^ + .^' 



29 8 144 

which is very close to 

« = cosP=i-g+^-..., 

as usually assumed. 

Inserting the v^ue of n found in the expression for M we get 



yI/=2|i-cos/8-2«^ry (cot ipyA 



which is just half the value of /* in the zenithal system. Thus the proper 
choice of the constant n has enabled us to find a projection which has only 
half the average error of a similar projection in that system. Gur map is 
however now fan-shaped instead of circular. 

Minimum Error Conical (approximately) Equal 
Area Projection. 

Taking now the second part of our general solution we have 

r = 2 (cot I6f f(tan ^eyde + 2B (cot ^Of, 

and on forming the expression for M and differentiating for the value of B 
to make this a minimum we shall find we must take 

^= {(tan |a)» ^^ - (tan i;8)» /J}/{(tan W - (tan W) 
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as before in the general case, and on substituting this we shall find as before 
that M reduces to 

OS a - cos ;8 - 2« C {^^ I {{tsxi ^/?)» - (tan Ja)^}] . 

The value of n will be the same as in the Orthomorphic system whether 
the map extends from the centre or over a zone a to /3, i.e. 

B^ & 
8 144 

m the first case and = cos v I i h . . . ) in the second. 

\ 30 / 

And when we form the expression for M in the general case 

r = A (tan ^df + (cot ^Of ( f(tan \efd0 + £\ , 

it will be found that there are no products of the several terms, and 
therefore the same values of n will apply ; also the value of M will be 

half the value just found, viz., M = -^-—? for map extending to centre and 

, _ 28^ sin x . 

M= for a zone. 

45 
It therefore appears that the Minimum Error Conical Projection is 
simply the arithmetical mean of the Minimum Error Orthomorphic and 
(approximately) Equal Area Projections and that all three give a map of 
the same total size and limiting radii, the intermediate radii in the two 
constituents being arranged to give orthomorphisrh oij, equality of area as 
the case may be and the absolute minimum error radii being the mean of 
the two. 



Numerical example of Minimum Error Map of a Zone. 

We will now illustrate our results by computing the constants and scale 
errors for a map extending from a = 30° to /3 = 90° from the pole of 
projection. In this case x = 6o° and cosx = i, and we shall adopt this 

value for n instead of the more correct value cos x ( ^ "^ •■• ) f^"" 

simplicity, and also because the term in 8^ does not affect the first term 
in the expansion for M. 

It will first be necessary to evaluate the function I (tan \BY dO and this 

we shall do, so as to give seven equidistant values between 30° and 90°, by 

Weddle's Rule. This Rule for evaluating the definite integral I. u^dx 

from seven equidistant ordinates u^, u^, u^, ... Ug is 

3^ 



r 



u^dx = — {ii„ + u^-^- Ui + Ui; + U2 + ^ {ui + U3 + Us)\> 
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where h is the value of the interval between abscissae. It is absolutely 
accurate if 5 th differences are constant, and errs in excess by only y^^A'^,, 

if 6th differences are constant. Applying this to evaluate I (tan ^6)^d6 for 
intervals of 10° we tabulate tan ^0 for intervals of 0° 50' as under : 





e 


i« 


tanj9 


« = (tanp)* 







0° 0' 


0° 0' 


O'OOOO 


o-cx)oo 




I 


i°4o' 


o°5o' 


■0145 




•1204 


2 


3° 20' 


i°4o' 


■0291 


•1706 




^ 


5° 0' 


2° 30' 


•0437 


•2090 


■2090 


4 


6° 40' 


3° 20' 


•0582 


■2412 




■; 


8° 20' 


4° 10' 


•0728 




•2698 


6 


10° 0' 


5° 0' 


•0875 


•2958 





Sum -9166 -5992 

5 X '5992 2-9960 

3-9126 

— X 3-9126=1-1738. The interval /z = i°4o' = cm. 

10 -^^ '-^ ^ 34377 

, 1-1738x100 

Integral= — — =0-0341. 

^ 34377 

Applying this method to successive intervals of 10° we get the following 
values of the integral. 



Interval 


|(tanjfl)*rfi9 




Total from o' 


0° — 10° 


0-0341 


10° 


0-0341 


10° — 20° 


-0630 


20° 


-0971 


20° — 30° 


•0821 


30° 


-1792 


30°— 40° 


■0980 


40° 


•2772 


40°— 5°° 


-1123 


50° 


•389s 


50° — 60° 


•1259 


60° 


■5154 


60° — 70° 


■1393 


70° 


•6547 


70° — 80° 


•1529 


80° 


-8076 


80° — 90° 


•1671 


. 90° 


■9747 


go° — 100° 


•1824 


100° 


1-1571 


100° — 110° 


•1994 


110° 


1-3565 


110° — 120' 


•2188 


120° 


i'5753 



re 
There is one case when the definite integral I (tan ^0)''d6 can be expressed 

in known functions and that is when 6 = -, when it can be expressed by 

2 

the 1/' functions of Gauss as - lij/ /^^— ^j - i/r (- — -ji . In this case since 

n = ^it becomes | {i/' (I) - i/' (|)} = 0-9749 as against 0-9747 in our Table. 
As this is only computed to 4 figures and depends on several additions 
the agreement is close enough. Using these results we shall find ^ = i -0870, 
.5= 0-1129, and we can tabulate our results as follows : 
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The first two columns are obvious. 
The third contains the values of Itan \BdB just found and underneath 

is written the sum I + ^, as this is to be multiplied by (cot \&f' as shewn 

in the fourth column. 

Columns 5, 6 and 7 are obvious. Column 8 is obtained by taking the 
differences of columns 4 and 5 and multiplying the result by nrjsin 0^ since 

^''/^^ - ' = ^ {^ ('^" *^)" - (^°' i^)" ( r ^ ^)} • 

These differences are set under the figures in column 4. 

Column 9 contains the squares of the scale errors, those of the radial 
scale being set under those of the parallel scale. 

Columns 10 and n are obvious. 

Summing column 1 1 by Weddle's Rule we find 



- _ o'oo8q67 X 10 

M= ^— ^ = o-ooi [;6. 

57-3° 



28= 



The formula Af= '^ gives 0-00151, the difference in the last figure 

bemg probably due to taking n = cos x instead of cos x ( i h ) ■ Dividing 

by the area of the zone we find /a = o"o42. Thus we have a projection in 
which the maximum scale error is + i2'5 per cent, on the upper parallel, 
+ 8-7 per cent, on the lower parallel and — 4'S per cent, on the mid-parallel, 
*hile the mean square scale error is 4-2 per cent. The standard parallels 
will be near 40° and 80°. 

Map of Hemisphere in Minimum Error Conical System. 

We will now calculate the particulars of the map of a hemisphere. 
Since a = o we put £ = and taking n = cos J/3 = ■o-7o7 1 we find A = 0-8369. 
We calculate the following Table (p. 30) in which the interval has been 
taken as 15°, in order to apply Weddle's Rule. 

We have thus a projection of the hemisphere in which the maximum 
and mean square scale errors are just about half those in Airy's Zenithal 
system. There is one standard parallel at about 65° from the pole. The 
map will be fan-shaped with a gap = 27r (i - n) = 104° about. This shape of 
map might not be so convenient for some purposes as a circular map of a 
hemisphere. 

Simplification of formulae by making radial scale true. 

An inspection of the scale errors in the Tables will shew that those 
perpendicular to greatly exceed those along the radii, particularly in the 
case of the zone map, where the radial scale is true at each boundary. We 
shall therefore simplify our formulae and sacrifice very little in accuracy if 
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we assume that the radial scale is true always, and then determine our 
constants so as to make the perpendicular scale error squared a minimum. 
We write therefore r = a + Q, where a is a constant, then 
dr\dQ = I, nr I sin 6 = n{a + $) sin 6, 

M= \ {nr/sin $- if sin 6d6= j cos a - cos |8 - 2« (y8 - a) (a + ^ (/3 + a)) 

„/,, tan^^S f^6d6 rP6''dO\] 

-vn^id" log ; — ^ + 2a -^— + ^— . ) I . 
V ° tan ^a ja, sm Q J^ sm 0j) 

Unfortunately the integrals 16 cosec 0d6 and 1 6^ cosec 6 dO cannot be ex- 
pressed in known functions. They could if necessary be evaluated by their 
expansions or by mechanical integration, but calling them for a moment 
Si and ^2 and putting dMjdn = o, we find the turning value 

a + hiP-"-) 



i = (/3-a). 



„, tan A/3 o c- 

tan ^a 



Putting next dMjda = o, we find 

13- a- nSi 



, /tan i/3\ ' 

Substituting this in the expression for n, we find after some reduction 

_ ^ (^ + g) (log tan ^jS - log tan |a) - S^ 

^''^ "''' ^^(logtanJyS-logtania)-^!^ ' 

whence M = cos a — cos (3 — n {j3 — a) (a + ^ {/S + a)). 

If the map extends to the centre then a = o and a = o. The quantity 

log tan 1/8 /log tan I a becomes infinite, so dividing the numerator and 

denominator of the right hand side of above expression for n by it, we find 

n=^^yS,; M=i-cosp-iP'IS^. 

Now S,= 6' cosec e d6 = '--+ '^ ■ '^ ' 

Jo 



2 24 2160 



Jo 
therefore 

n=i-i^+ -4 — +..., 

12 2160 90720 

J/_^' ^\ ^ ^" I ... ^'1 ^' ^' I -^ 

2 24 720 40320 '" 2 24 4320 181,440 

^^/_^ \ 

864 V 60'" J' 

which is 2/5ths the value of M in the Equidistant Zenithal system. 

If the map extends over a zone from ^ = a to = 13, the expressions for 

n, « and M can be simplified by expanding as follows, for which method 

the writer is indebted to Dr T. J. I'A. Bromwich, F.R.S. Putting 

^{P + a)=x and i{/3-a) = S, 
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e 



let i/sin (x - e) = A + ^i£ + ^2 — + ^3— + ..., 

so that 

>/ - ^_ . ^ _- t^os X . . _ I + cos^ X . cosx(5+c°s^x) . 

^"-sitix' ^'~^^x' "~ sin^X ' ' ^i^^'x ' 

and so on; further write 

Z= -T—. r=2S \A^+ A^ — + Ai + ... , 

;_{ sin(x-e) V 6 I20 /' 

j_s sin(x-£) 3 \ lo 280 ) 

,^ /■+« cVe 2SV ^ ^ 8' ^ SO \ 

j-a sin(x-€) 3 V -^ '10 '56 ;> 

28i}/ „ , . . J , . 

then « = ytt — 7^2 • Substituting and reducing we find 

«=cosx (^i +— ...j, a = tanx(i --)-x. 

(^ 8' S= „ 8^ 8=\ . S= 

7>/=2sinx 8-— + ...— S+-r+ — = 2sinY — , 

'^ \ 6 120 6 72/ '^ 45 ' 

which values for n and M agree with those found for the Minimum Error 
solution so far as our expansions go. We are therefore justified in rejecting 
the rather complex formulae of that solution and adopting the uniform 
radial scale. 

We will now illustrate our results by calculating the particulars of the 
same zonal map as before, extending from 30° to 90° from the pole of 
projection. Putting 8 = 30° = 0-5236 cm. and x = 6o°j we find = 0-5899, 
« = 0-5046, and our Table stands as under : 



I 


2 


3 


4 


5 


6 


7 


8 


6 


B in cm. 
9 + a 


sine 


nrjsin 8 


nrjsin 8-1 


Squares 


VsinS 


nrjsm 8-1 

in Minimum 

Error 

System 


30° 


0-5236 
I'li35 


0-5000 


I-1237 


+ -1237 


0-015301 


0-007650 


+ 0-1252 


40° 


-6981 
1-2880 


-6428 


roiio 


+ -0110 


121 


78 


+ -0130 


50° 


-8727 
1-4626 


-7660 


"9634 


- -0366 


1340 


1026 


- -0364 


60° 


1-0472 
1-6371 


-8660 


-9540 


- -0460 


2II6 


1832 


- -0452 


70° 


I -221 7 
1-8116 


•9397 


•9728 


- -0272 


740 


69s 


- -0280 


80° 


13963 
1-9862 


-9848 


1-0178 


+ •0178 


3«7 


312 


+ -01 60 


90° 


♦ 1-5708 
2-1607 


I -oooo 


1-0912 


+ -0902 


8136 


8136 


+ -0870 
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The columns are self-explanatory and in the last are given the perpen- 
dicular scale errors in the absolute Minimum Error Projection, for com- 
parison with those in column 5. By Weddle's Rule we find J/'=o'ooi59 
as against o"ooi56 in that system. The formula to first term gives 
M=o'oo\^\. We see that the map has two standard parallels and that the 
scale' error at the mid-parallel 60° is 

«(a + x)/smx-i=(i+-j(^i--j=-^. 

To find approximately where the standard parallels lie, let y be their angle 
from the centre, then « (a -)- x — y) = sin (x — 7), therefore 

cosx(i "^^)('^^"x(^ ~l) -v) = sinx(i-iy')-ycosx, 

whence 7 = ± -^S = + o'588, 

v3 

that is, the standard parallels lie at approximately one-fifth the depth of the 
map from the upper and lower edges. 

It will be observed that the mean square scale error is independent of 
the distance from the pole of projection, or what we may call the mean 
co-latitude of the map, using co-latitude in a general sense as angle from 
the pole. This is an important result, as it affords us great liberty of choice 
in fitting any particular country to a map in this system. We have only to 
choose our pole so that the depth of the map 28 shall be as small as 
possible, and as the choice of mean co-latitude x enables us to select 
almost any curvature for our parallels, we can suit almost any shape of 
country. Hammer has published Tables to facilitate the construction of 
Conical Projections with pole in any position, but when the pole has a 
skew-position it will be found that generally the meridians and parallels are 
lines of contrary flexure and the map has an ugly appearance. Practically 
speaking the only useful positions are when the pole corresponds with the 
terrestrial poles or when it is at an infinite distance in the perpendicular 
direction, leading to transverse cylindrical projections. 



Minimum Error Plate-Carree Projection. 

When the mid-parallel in the Equidistant Conical Projection is of co- 
latitude 90° the projection becomes cylindrical, and if the pole of projection 
is at infinity in the direction of the terrestrial poles it is called the Plate^ 
Carree Projection. When, however, the pole is at infinity in the transverse 
direction the projection becomes the Cassini Projection. In both cases can 
these projections be improved by use of a suitable scale constant. 

Considering the Plate-Carree Projection, let us now count the angle 6 
from the equator instead of from the pole. The quantity n becomes zero. 
The parallel scale error, if the cylinder is tangent at the equator, is repre- 

Y. 3 
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sented by sec 5- i. Let now A be sl scale constant and let us determine 
A so that 

2 ( (A seed- lY cos OdO = 2 f— log ^"^^!"! - 2^8 + sin sl 

shall be a minimum. We shall find 

A = 2S/{log (1 + sin 8) - log (i - sin 8)} = I — -^ ... ; 

28^* 8^ 

J/= 2 [sin8- 287{log (i +sin8) — log (i — sin 8)}] = — and /"■ = — 7 — ■ 

These are precisely the results which we should have obtained by making 
X = 90° in our general formula. 

If we take A= i, or make the cylinder tangent at the equator, 

81 8' 8^ 



,. r • c. c 1 1 I +sm 8) 8' 

J/ = 2 -^smS— 26 + i log -. — sJ- = — ; «,= 

( '^ ° i-sm 8J 10' "^ 



V20' 



By creating two standard parallels instead of one, we have introduced a 
negative scale error along the middle parallel of — 7- , we have reduced the 

mean square scale error in the ratio ^45 : ^20 = 67 : 4'5, and the maximum 

8^ 
scale at the extreme parallels has been reduced from sec 8 = i + — . . . to 



?)(-!>-: 



Minimum Error Cassini Projection. 

Tissot has pointed out (p. 162) that the Cassini Projection is simply the 

Plate-Carree turned through a right angle, so that what was the equator 

becomes the principal meridian. Our result above shews that by reducing 

the scale along this meridian by one-third the scale error at the extreme 

east or west edge of the map, creating what we may call two standard 

pseudo-meridians (or sub-parallels), we reduce the mean square error of the 

map to a minimum. A device of this sort is used in the projection of the 

International Map on the scale of i/million. The scale error on the 

principal meridian is made equal and opposite to that on the extremes, i.e. 

8- 8^ 

the scale constant used is i instead of i - — . 

4 6 



Minimum Error Mercator Projection. 

The work we have just done enables us to find the Minimum Error 
Mercator Projection which is simply the Cylindrical Orthomorphic. In this 
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the radial or meridian scale instead of being uniform is the same as the 

parallel scale- Using (r — rw) as distance from equator, we have 

2 

d{r- r^)/d^ = sec e, 



which gives (r — r„) = log tan (- + AS) = i log -. — ;; . 

The expression for M is simply double that in the Plate-Carrie, 
A is the same, but must be applied to the perpendicular as 

meridian coordinates, and fi. is J2 times as great = ./ — • S^ 



Minimum Error Conformal Projection of Gauss. 

As the Cassini is simply the Plate-Carree turned through a right angle, 
so the Conformal Projection of Gauss, used for the survey of Hanover and 
more recently for the survey of Egypt, is simply the Mercator applied to a 
meridian. We see therefore that its mean square scale error can be reduced, 

by multiplying both coordinates by a factor A = i ..., in the ratio 

6 72 

671 : 474 and its maximum scale error in the ratio 3:2. The limiting 

depth of the rectified Mercator or width of the rectified Gauss Conformal 

is ^ ^ log : — ^ = S, the same as in the Plate-Carree or Cassini-Soldner. 



I 



■ sin 8 



We shall now compare our Minimum Error Equidistant Conical Pro- 
jection with some others of the same class. The first of these is : 

Euler's Map for Russia with two standard parallels. 

This projection is described in detail in Germain, pp. 179 — 183. It has 
the property that the absolute errors on the extreme parallels are equal, and 
also equal to the error of opposite sign on the mid-parallel. This leads in 
our notation to the two equations 

^ (« + X — ^) ~ sin (x — 8) = sin X - « (« + x) = '^ ('^ + X + ^) ~ sin (x + 8), 
whence n = cos x sin 8/8 = cos x ( i ~ t + ■ • • ) i 

a + X = tanxiScot|8 = tanx(i " i^"-^---)- 

The map extends from co-latitude 20° to co-latitude 50° and Euler finds 
by a calculation not so direct as ours « =0-8098, « = o'o853 = 4° 53' 22". 
As this value of a is so nearly 5° = o'o873 Euler adopts the latter value and 
so fails to reproduce one of his conditions, viz., that of equality of errors at 
middle and extremes. The following Tables give the particulars of this 
map calculated from the Minimum Error formula, from Euler's own con- 
stants, and from the correct values of Euler's constants. 

3—2 
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Minimum Error Formulae. 





« = cos X ( 


I + — =o' 
30/ 


82156, 


a = tan x ( 


,.--^x= 


: 0-07927. 


e 


6 in cm. and r sin 9 


nrjsm 6 


nrjsm 8-1 


V= Squares 


VAnB 


20° 
25° 
30° 


0^34907 
■42834 

■43633 

•51560 

•52360 

•60287 

•61086 

•69013 . 

•69813 

■77740 

•78540 

•86467 

•87266 

■95193 


©•34202 
•42262 
•50000 


ro286 

I ^002 2 

•9906 


+ •0286 
+ ^002 2 

- ^0094 


o^oooSiS 

4 
88 


0-000280 
2 

44 


35° 
40° 

45° 


■57358 
•64279 
•7071 1 


•9884 

■9936 

roo45 


-•01 16 
-•0064 

+■0045 


134 
41 
20 


77 
26 

14 


50° 


•76604 


I -0208 


+ ^0208 


433 


332 



A/'= 0-0000320 by Weddle's Rule 
=0-0000313 by formula. 

Same Map with Euler's constants, ;^ = 0-8098, a = o-o8j27. 

o^ooio3o 352 + (o^oiio) 



20 


©•34907 
■43634 


0'34202 


I -032 1 


+-032I 


25° 


■43633 
•52360 


-42262 


1-0032 


+ -0032 


30° 


•52360 
•61087 


-50000 


■9895 


- -0105 


35° 


•61086 
•69813 


■57358 


•9857 


- -0143 


40° 


■69813 
•78540 


-64279 


•9894 


- -0106 


45° 


•78540 
•87267 


-707 1 1 


■9994 


--0006 


50° 


•87267 
■95993 


-76604 


I -0147 


+-OI47 



10 


4 


no 


55 


204 


117 -(0-0082) 


112 


72 








216 


165 +(0-0112) 



Af= 0-00003 5 2 by Weddle's Rule. 

If Euler's assumption had been reproduced we should have the scale errors 
at extreme and mid- parallels, multiplied by respective sines of 6, equal, 
whereas we have the figures shewn in brackets in the last column. 

Same Map with Euler's correct constants, « = 0-8098, a = 0-08535. 
8 Sin cm. and r sinff nrjsin8 nrjsmd-i F= Squares Vsmd 

10° 0-34907 0-34202 1-0285 +0-0285 0-000812 0-000278 + (0-00976) 





■43442 












25° 


■43634 
-52169 


■42262 


■9996 


- -0004 








30° 


-52360 
-60895 


-50000 


-9863 


- ^0137 


185 


93 


35° 


-61086 
-69621 


■57358 


-9830 


- -0170 


289 


166 -(0-00975) 


40° 


-69813 
■78348 


-64279 


•9871 


- -0129 


166 


107 


45° 


-78540 
-87075 


■707 1 1 


-9972 


- -0028 


8 


6 


50° 


-87266 
-95801 


-76604 


1-0128 


+ -0128 


164 


1 26 + (0-00980) 



Af=o-oooo427 by Weddle's Rule 
= 0-000041 by formula. 
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The assumptions made are how very nearly reproduced, but, as we see, 
Euler by using erroneous values of his constants, was fortunate enough to 
reduce his mean scale error. 

To calculate the formula for Tkf or /«. we should strictly speaking substi- 
tute the values of n and a in the equation 

J/ = cosa-cos/3 + «^( a'"log : — f- + 2a 6cosec6dd+ I ff'cosec 6d0] 
\ ° tan ^a Jo. A / 

-2«(/3-a)(a + |(j8 + a)), 
but as we know /a is to the first term independent of the co-latitude we are 

TT 

entitled to take the value for y = - and write 

2 

M=2 jsin 8-2^8 + ^^4° log _ . A sinx 

Now A we know is the scale at mid-parallel 

^ 8* 
= n(a + x)/sm x = J sin 8 cot ^8 = i - - + -^ . 

Substituting this we find 



120 '^' 1^ /^ \240/ 



Hinks, pp. 83 — 87 (following Close, Textbook of Topographical Survey- 
ing, p. 107 et seq.), considers five different assumptions for this projection 
including Euler's (Case I). We will now exemplify his Case III in which 
the Scale Errors (not the absolute errors) along extreme parallels are made 
equal and opposite to that along the mid-parallel. The maximum negative 
scale error does not occur exactly on the mid-parallel, except when x = 9°°, 
but it is very near to it, so that Hinks' Case V is practically the same as 
his Case III. We have then the following two equations 

i^ (« + X - 8) _\-n{a + ^ _n{a-\-x + ^) 
sin(x-S) sinx sin(x + S) ' 

whence we get 

« = cosx- 2 tan 18/8 = cosx(^i +7^+7^-"j' 

/ 8^ 8* \ 
a -h X = tan x • 8 cot 8 = tan x(i ■••)• 

8^ 8^ 
Scale at mid-parallel = 2 tan i8 cot 8 = i ... the same as m Euler's 

map ; the value of M in this projection will therefore be the same as in 

Euler's to first term = ^ ^^^X -^yg ^^^ n = o-?,2i?,; a = 0-0733, and Table 
120 

as follows : 
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9 in cm. and r 


sine 


nrjsm 


nrlsin 0~ i 


V= Squares 


V5\n0 


20° 


0-3491 
•4224 


0-3420 


I-0174 


+ 0-0174 


0-000303 


0-000104 


25° 


■4363 
■5096 

■5236 
•5969 


•4226 


•9934 


- -0066 


43 


18 


30° 


•5000 


•9834 


— -0166 


275 


137 


35° 


•6 log 
•6842 


■5736 


-9826 


- -0174 


303 


174 


40° 


•6981 
7714 


-6428 


-9886 


— -01 14 


130 


84 


45° 


7854 
•8587 


-7071 


1-0003 


+ -0003 








50° 


■8727 
•9460 


•7660 


roi73 


+ -0173 


299 


229 



7J/= 0-000044 by Weddle's Rule 
= 0-000041 by formula. 

The position of the standard parallels in this, and in all cases where 
the scale error on mid-parallel is equal and opposite to those on extreme 
parallels, is given by 0^ = i^^ ^•^- about one-seventh the depth of the map 
from edges. 

We will now discuss two remarkable projections of this class invented 
by the Rev. Patrick Murdoch, F.R.S., and described by him in a paper 
published in the Philosophical Transactions of the Royal Society for 1758. 
But first we will investigate the condition that the total area of the map 
shall be equal to the area on the sphere which it represents. The area of a 
zone on the sphere between co-latitudes a and ^ is 

2tTB^ (cos a - cos ^) = 27r^^ 2 sin 8 sin x- 
The corresponding area on the map is 

rnr {r^ — rj) = mr {r^ + r^ {r^ — r^) = irnrE? . 28 (a -t- x)- 
Therefore for equality of areas «8 (a->ry) = sin 8 sin x, or 

8^ h* 
n (a + y)ls'm y = sin 8/8 = i — ;r H .... 

^ ^'1 ^ I 5 j;20 

Thus the Minimum Error Conical Projection has to a first approximation 
its total area true, as in the zenithal system. 

Murdoch's First Projection. 

Murdoch's Paper is very concise and he does not give the full proof of 
his result. He takes n = cos x simply, and then gives a rule for drawing the 
map which in effect amounts to taking a + x = tan x sin 8/S, and this we see 
is the necessary condition. 

Murdoch's result was afterwards taken up by the German mathematicians, 
Tobias Mayer and Albers (see Germain, p. 185). The fornier professed to 
have discovered some flaw in Murdoch's argument, and the latter elaborated 
the result and also developed his own Equal Area Projection of this class 
in which the radial scale is of course no longer uniform. The writer has 
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19 



examined Murdoch's original Paper and can discover no flaw in his state- 
ments j and it is remarkable that Murdoch by a kind of intuition hit upon 
a projection which has to a first approximation the same mean square scale 
error as the absolute Minimum Error Projection of its class. 

• -7 H m the formula 



For if we substitute A= i 



f S^ 



2 A + iA^ 



S'M 



8^ 
V 6 24 



smx, 



we find M= — sin v as before. 
45 ^ 
To exemplify this we will calculate the map from co-latitude 30° to 
co-latitude 90° on Murdoch's First Projection. We shall find 
n = cos 60° = o'5ooo, a = tan x sin S/8 — x = o'6o68. 



e 


8 in cm. and r 


sinfl 


nrlsin 


nrjsin B - i 


F= Squares 


Vsm.e 


30° 


0-5236 


0-5000 


I -1304 


-ho- 1 304 


0-017004 


0-008502 




1-1304 














40" 


■6981 
'■3049 


•6428 


I -01 50 


+ 


•0150 


225 


145 


50° 


•8727 
1-4795 


-7660 


•9658 


— 


•0342 


1 170 


896 


60° 


1-0472 
1-6540 


-8660 


-9550 


-' 


-0450 


2025 


1754 


70° 


1-2217 
1-8285 


■9397 


•972S 


— 


-0372 


1384 


1300 


80° 


1-3963 
2-0031 


•9848 


10170 


-f 


•0170 


289 


285 


90° 


1-5708 
2-1776 


i-oooo 


1-0888 


-1- 


•0888 


7885 


7885 






M= 


0-00164 by Weddle's Rule 
0-00151 by formula. 







We thus find scale errors practically identical with those of the absolute 
Minimum Error Projection of this class, and a mean square scale error 
greater in the ratio ^164 : ^^156 only, and this for a zone map of a very 
large extent. The following rather disparaging remark of Germain (p. 186) 
is therefore not justified : " Quelle que soit la m^thode adoptee, cette 
projection ne pent convenir pour une bande de plus de 8 ou 10 degres de 
latitude, et dans des limites aussi rapprochees presque toutes les projections 
coniques ont a peu pres les memes avantages ; quant a la question de 
conserver la superficie totale de la carte, elle n'est que secondaire et oblige 
au contraire a comprimer les superficies et les distances dans la zone 
moyenne, tant en longueur qu'en largeur, et a les etendre dans les deux 
zones extremes dans toutes les directions." The statement as to the com- 
pression and extension of the depth of the map is certainly not correct, as 
the radial scale, is everywhere uniform and true. 

Murdoch's Paper was followed a month later by a letter from Mr Wm. 
Mountaine, F.R.S., on "A short Dissertation on Maps and Charts," in 
which he challenges Murdoch " to draw a map in which the superficies of 



4° 
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a given zone shall be equal to the zone on the sphere, while at the same 
time the projection from the centre is strictly geometrical." 

Murdoch replies in a Note printed as an Addendum to Mountaine's 
Paper giving the Rule for his Second Projection which is a Perspective 
Conical Projection, and he says, " But here the degrees of latitude towards 
the middle will fall short of their just quantity and at the extremities exceed 
it : which hurts the eye." He then indicates the Rule for 



Murdoch's Third Projection. 

In this he takes a + x = tan x ■ ^ cot 8. This value of a + x is quoted by 
both Germain and Tissot, but it does not appear to have been noticed that 
the condition of equality of total area exacts a corresponding change in the 
value of n, viz., 



sin S sin X _ cos x sin 8 



a+x 



a + x-- 



( 8^ S^\ 

tan X I • 

V 3 45/ 



8^ cot 8 
45^ 



cosx 



6 3 
— + - 
6 36 






Scale at mid-parallel = i - ^r h , 

6 120 



; M= — sin x as before. Applying 
45 



these results to the same map as before we get « = o'5265, (3; = o'645i, and 
our Table stands as under : 



e 


B in cm. and r 


sin« 


nrjs,m 8 nrj^m d - 1 


V= squares 


Ksine 


30° 


0-5236 


0-5000 


I-I027 +0-1027 


0-010547 


-005273 




1-0472 












40° 


-6981 
I -22 1 7 


•6428 


I -0006 + 


-0006 








5°° 


-8727 
1-3963 


-7660 


-9596 


•0404 


1632 


1250 


60° 


1-0472 
1-5708 


•8660 


-9550 


-0450 


2025 


1754 


70° 


1-2217 
17453 


•9397 


-9778 


-0222 


493 


463 


80° 


1-3963 
1-9199 


-9848 


I -0264 4- 


-0264 


697 


686 


90° 


1-5708 
2-0944 


I -0000 


rio27 + 


•1027 


10547 


10547 






M= 


0-00165 by Weddle's Rule 










= 


0-0015 1 by formu 


la. 







We thus have a map which has the following properties : 
(i) The total area is true, 
(ii) The radial scale is everywhere true. 

(iii) Two parallels are true to scale, distant about |^th depth from edges, 
(iv) The scale errors at extreme parallels are equal and about double 

the scale error of opposite sign at mid-parallel, 
(v) The average scale error is to first approximation the same as in the 

Minimum Error Projection of this class. 
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{vi) The maximum scale error is less by Jth than the same quantity in 

the latter. 
By substituting the values of n and a we find 

Scale at extreme parallels = tan 8/8 = i h — h . . . , 

32 8* 

,, mid-parallel =sin8/S=i- — h 

6 120 

This Third Projection of Murdoch therefore would appear to be on the 
■whole the very best of the Conical Class, better for all practical purposes 
than the absolute Minimum Error, as the reduction in the maximum error 
more than compensates for the slight increase in average error. It is very 
easily computed whether by the full formulae or by the expansions in terms 
of 81 This is a curious reflection on Germain's remark, p. 189, " Ce 
troisieme systeme etant, sous tous les rapports, inferieur aux deux premiers, 
nous ne nous y arreterons pas." It is also some satisfaction that it was 
invented by a fellow-countryman, who though born and educated in Scotland, 
was Rector of Stradishall, Suffolk, England. 



Equidistant Conical Projection with total area true. 

As in the Zenithal System the condition that the total area of the map 
shall be true introduces practically no difference. It affects the value of 
-Af only in the third term of its expansion, and for n we have the equation 

27r I nQM = 277 (r — cos p), 
h 
whence 

« = 2 (i -cos^)//8^= 4sin==i/3/^^= I --^ + J^- ..., 

■differing in the third term only. 

Substituting this in the expression for M we find 

M=i-cosS-nBr' + n^ A" coi&c 6 d6 = -^ ( i - •=- ... 
Jo 864 \ 60 

as before. 

The maximum scale error at the boundary 

= «/3/sin/3 = 2 tan J/8/^= i + ^ + 7^ - 
as against 

\ 12 2160/ \ 5 360/ 12 2100 

The following Table exhibits the map of the hemisphere in this system, in 
which n = 0797 : 



42 MINIMUM ERROR 



e 


r 


sin 9 


nrjsin 8 


nrjAa 8-1 V 


= Squares 


Fsin8 


o° 


o-oooo 


o-oooo 


0-7970 


- -2030 1 


D-04I2I 


O'OOOOO 


15° 


•2618 


-2588 


•8062 


--1938 


-03756 


972 


30° 


■5236 


-5000 


•8346 


--1654 


-02736 


1368 


45° 


7854 


-7071 


■8852 


-•1148 


■01318 


932 


60° 


1-0472 


-8660 


■9639 


— -0361 


■00130 


113 


75° 


1-3090 


■9659 


I -0801 


+ •0801 


■00642 


620 


90° 


1-5708 


I -0000 


1-2519 


+ -2519 


•06345 


6345 






M= 


0-0168 by Weddle's Rule 










= 


0-0167 by formula. 







In the Minimum Error Projection of this class the value of M was 
0-0091 and the value of the scale error at /3 = go° was +o'i835, so the 
merits of the uniform radial scale are not so near those of the absolute 
minimum, in so large a map as a hemisphere, in this class as in the 
Zenithal. It is still, however, much superior to Airy's for the hemisphere, 
in which Afyfus 0-0391 and the scale error at the boundary was 0-3863. 

It may be objected to this fan-shaped Conical system that it makes 
a rather ugly map when n is grisater than say ^, as there is a wedge-shaped 
gap = (i — ?«) 360° not filled up. This gap may, however, be filled up by 
continuing the map, which has the advantage of placing in juxtaposition 
those parts which abut on the boundaries, so that continuous measurements 
may be taken across those regions. Theoretically speaking this advantage 
holds however small the limiting radius /8 of the map, but it is doubtful if 
there would be any advantage when the gap is less than 30°, or p less than 
about 60°, for all bearings from the centre are too small by about one- 
twelfth and, when the projection is so near a complete circle, this might 
be inconvenient. 

Plate-Carree and Mercator Projections with total area true. 

In these cases we have Ap = I cos Odd = sin /8 and therefore 

Jo 

B^ B* 
^ = sin/3/fi=i-^4--^... 
6 120 

as against ■■■ ^ V ■ ■ - 

in the absolute minimum system. It will be found that the value of Mia- 
not altered in the first two terms of its expansion. 

Expansion of General Formulae in case of Zone Map. 

In the case of the Zone Map the limiting radii are the same, whether 
the map be Orthomorphic, approximately Equal Area, or the mean of the 
two, viz., 

(tan ^/S)" ( (tan ^efdO 2 tan ia ( (tan 16)" dO 

Ja Ja 



2 



(tan i/3)« - (tan ^af ' " (tan ^y8)»- (tan ^af ' 
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Bep-h Of ™.p= ,/;(^ 1.,.. i <.rff:;;::;|:p =.-.. 

Now we have found 

f (tan je)ve = .8 (tan Jx)" (i + ? ^J^tZ^lA 
la. ^ " ' \ b sin'^x 

8* n*' — dti? cos X + 7«' cos^ x + 4^° ~ 5^ cos x — ^ cos^ x 1 

120 sin*x J 

= 28 (tan Ix)" ] I -f- cot^ X f if 'z = cos x ( i + — ) 

and =28 (tan |x)" ] i cot^ x f if ^^ = cos x simply. 

We have also found that 

(tan m- - (tan ^a)- = ^^'^^l^^^xr f ^ 1 4^-- 6.cos x+ i +cos'x 
^ 'i'^i ^ "i I sinx 1 6 sm=x 

i6«*— 8o;^'cosx + 4°'2''+ loo'z^cos^X" To'^cosx — 3o«cos^X 1 
8^ + i8cos^X + cos*x+ 5 h. 



1 20 sin^ X 1 

and replacing 2» by « we find 

(tanW-(tan^ar= ^"^^^""^^^7 i+g ^'-3""°^/,^^^"°^'>^ 
^ ^'^' ^ ^ '^ sinx 1 6 sin^x 

n* — io;z'cosx+ io«^+ 25«- cos^x "35^^ cosx— 15^2 cos^x 1 

8^ + i8cos^x + cos''x + 5 \ 

120 sin^X "J 

Dividing out we find 
{(tan i^)^ - (tan \ a)»'}/{(tan J^)" - (tan |a)»} 

_ (cot l^x)" / 2^ 3«^- 3^2 cosx 

2 I 6 sin^x 

8** 75«^— 9o;7^cosx-6o«^— i5«^cos^X+ 7S'^cosx + i5« cos'x 1 

360 sin^X > 

J_cotix)!L_Sl^„,. S^cot^XJ 
2 (.60 '^ 24 J 

(cotix)" r S'cot^xl -f '/ s'\ 

= -=^ ^^^^ \ I + ^f if « = cos X ( I + — 

2 I 40 J \ 30/ 

and 

(cotix)" f S^cot^x) -c -1 

= -^^ ^^=^ { I H ^}- if « = cos X simply. 

21 24 J 

if « = cos Y I H 

^ \ 30/ 



Therefore 

8^cot^X^ ■■■ / 8= 



1 -;-„ = 281 I + 



45 / ^ \ 3c 

8* cot^ 



./ 8*cot^x\ -r ■ , 

and =28(1 + j it n = cos x simply. 

These shew how nearly the depth of the map is true. 
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We shall now obtain an expression for the radius of the mid-parallel 
and it will be simpler to treat the orthdmorphic part separately, where 



r=2A (tan J6l)" = ["2 / ^ (tan i6»)''^6i/{(tan JyS)^'' - (tan ^af""] 
= (using above expansions) 



(taniS)" 



sinx (^ ^' 3;>2^-5^cosx+i+cos^x l /.„„ 1 m» 
Mtan^xr i 6 ^i^ |(tan^&) 

8^ /tanJ^'N" ., / 8' 



/ 8^ /tan|6'\'' .^ / S^\ 

= tan V I — ^ if « = cos Y I + — . 

^\ sAtanIx/ V 3°J' 

Thus when 6 = x, ^x = t^^" X ( i ) > 

the same as in the Equidistant System. 

Now put ^ = X + e, where e is the angle from mid-parallel to any 
other; tTien 

(tanie)"=(tani(x-He)r 

,^ n \„ r 'ze fr' — ncos\€' n" — vi'' cos -y + n + n cos' X ^ ) 

■ ■ =(tan|x)Mi + -^ — + ^2 ~- + -3 ^T+ ■■■( 

{ smx sm^x 2 sm'x ° J 

i-i-£C0tx+€C0tx — +^cotx...| if ;^ = cosx^i+— j, 
therefore 

This gives a very ready method of calculating and plotting the Minimum 
Error Orthomorphic Projection of a zone, when once the radius of the mid- 
parallel has been computed. 

In the approximate Equal Area Projection 

r=2 (cot i ey [(tan ^efdO + zB cot 1 6, 

' (cot^jS)^" f (ia.n ley dO- {cot laf"" r {ta.niefde 

where .S = ^5^-^ ^^ ■ 

(cot|a)^''-(coti/3)-^'' 

we shall find for the radius of mid-parallel 

"(cotiaf f^ + {cot ilSy'j^' 

r = 2 ("cot iv"!" ^ ^ 

'■x MCOt^x; L (cotla)^''-(cot^i3)^'' .■ 

Expanding this expression in the same way we shall find 

sinx / 8^ . / 8\ .„ / 8' 



''x = 



and evaluating ^e — f^ in the same way we get ^ [^ + t ~ t) l so that the 

mean of this and the orthomorphic part is simply r^ — r^ = e, the Equi- 
distant value. 
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The Albers Equal Area Projection. 

The Projection we have just considered is not strictly Equal Area since 
it satisfies the differential equation drjdd— i = — (nr /sin - i), whereas the 
strictly Equal Area Projection satisfies the differential equation 

dr/dd = sin 0/nr. 
The complete solution of the latter is 



•=( 



r„^ + ^sin^i(9] 



and is called the Albers Projection. When the Projection extends to the 
centre r=o when B = o, and r= 2sm.^6jjn simply. We shall now in- 
vestigate what value of n makes the mean square scale error a minimum. 

Putting Jn = I + «, so that — ^ = i — a + a^— . . ., we have 

r=(i~a + a' ...) (0 + ... ), 

^ ' V 24 1920 / 

^,./^e=,__ + __..._«(T-«)(^I-_...j, 

M= r \{drjde - if + (/-/sin 6 - i)^} sin OdO 

=n£(-3'^f'-"'(-S(-f)--'— '(-?)" 

whence for a minimum we shall find 

16 512 

or i+a=i 7+ , 

16 512' 

18^ B* 
and (,+^). = ^=i_^ + _^g, 

fi^ B^ 
which value of n differs very little from cos \P= 1 — ~ + -^ . . . proposed 

by M. Tissot to make the maximum scale error as small as possible. The 

/ /3^ B* \ 
limiting radius rp = /3li+'—^ fifi~-") ^"*^ '•^^ maximum scale error 

B^ 1 1 jS'' 
thereat = nr^jsm /5 = i + -7 H ^ which differ by small amounts in the 

third terms found for the approximate Equal Area Minimum Error system. 

The value for Mis -~ to the first term, the same as in the latter system. 
760 
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As to the Albers Projection for a zone, r = (r^^ + - sin^ \0] , 



the 



writer has not succeeded in deriving by rigorous mathematical processes 
the values of the constants r^ and n which render its mean square scale 
error a minimum, the requisite integrations and expansions being very 
complex. By writing the formula in the form r =^ {a {c- cos 6)Y> where 
a = 2ln and a {c— i) = r^, we have 

drjde = a sin 6/2 {a (c - cos 6)]K fir/sin 6 = 2 {a (c- cos (9)}V« sin 0, 
and 

M. r( ^^i"^^ _ yasin^g ^4{ c-cos6) _ 4{c-cos6)i _^ ^^.^^ ^^_ 
}a U{c-cose) (c-cos6)i asmd Ja j 

We can now expand the integrals in terms of the middle value 
X = i (/^ + ") and the limits 8 = i (j8 - a), and we may put in the expansions 

a = 2/n = 2/cos X (i + ^)i ^- cos X = sin^ X (i +/)/2 cos x, 
where e and /are corrections depending on 8^ and higher terms, the values 
of which are to be found so as to make M a minimum. The writer has 
not taken the expansions beyond 8^ and so far he finds the equation 
£+/+ 3-8^ = must be satisfied. 

By analogy with the approximate Equal Area function, we may however 
assume that if we make the scales at upper and lower limits equal and the 
total depth of the map true, we shall not be far from the absolute minimum 
mean square error conditions. These assumptions have been worked out 
by Fiorini (p. 269) and amount to taking 

I + cos a cos /8, 2 8^ cot 8 



whence C— cosx 



2 cos X cos S ' cos X sin 8 ' 

sin^x(i-l) 2(i-|) 

2 cos X cos X 



S2 22 

thus e = + -T, f= and our equation e +/+ \^ = o is satisfied. We also 

find 

/ 8^ / l\ 

« = cosx(^i + gj, r^ = tanx(^i--j, 

the same values as in Murdoch's Third Projection. 

To exemplify this result we shall compute the same map from a = 30° 
to ;8 = 90° (Upper table). We find « = 0-5265, a = 37988, C=i'i547. 

Tissot (pp. 222 — 6) has investigated the values of the constants of the 
Albers Projection to make the maximum scale error as small as possible. 
This involves making the scales at upper and lower parallels equal and also 
equal in value to the maximum scale of opposite sign which occurs near 
the mid-parallel. We shall now compute the same map on Tissot's assump- 
tion, for comparison (Lower table). In his notation r= ( r„^ + - sin^^fl) ^ 
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cos e Squares 

e C-cose a(C-cose) r nrlsine drjdd nrlsin e - 1 drjcie - 1 and y VsinB 

30° o'866o I "097 ro48 i'io4 o'9o6 +'io4 -'094 -oiio •0099 

•2887 88 

•0198 
40° 7660 1-477 1-215 "995 I '005 --005 +-005 -CXXX3 -oooo 
•3887 o 

-CXXX3 

50° -6428 1-945 1-395 '959 i'o43 -'041 +'043 -0016 -0026 

•5119 18 

-0034 
60° -5000 2-487 1-577 -959 1-043 --041 +-043 -0016 -0029 

•6547 18 

•0034 
70° -3420 3-087 1-757 -984 i-oi6 --016 +-oi6 -0003 -0006 

■8127 • 3 

-0006 
■80° -1736 3-727 1-931 1-032 -969 +-032 --031 -0010 -0020 
-98 1 1 10 

•0020 
■go° -oooo 4-386 2-094 I-I02 -907 +-102 — -093 -0104 -0190 

1-1547 86 

■0190 

By Weddle's Rule Ar= 0-00311 as against 0-00300 for double Murdoch's 
Third Projection and /t = oo6o, so we may be sure we are very close to the 
conditions for absolute minimum of M. 



-(i-cos«) 
d I -cos 9 and ^^ r tzrjsaxd drjdd Errors Squares and K Ksin 9 

30° 0-1340 0-5360 1-075 1-075 0-9^0 +-075 0-00562 0-00526 
1-1550 — -070 490 

0-01052 
40° -2340 '9360 1-247 '970 I '03 1 —-030 90 120 

1-5550 +-031 c6 

186 
50° -3571 1-4284 1-431 -934 1-071 - -066 436 720 

2-0474 +-071 504 

940 
■60° -5000 2-0000 1-619 -935 1-069 -'065 422 778 

2-6190 +-069 476 

898 
70° -6580 3-6320 1-803 -959 1-043 --041 168 332 

3-2510 +-043 185 

* 3S3 
80° -8264 3-3056 I -981 1-006 '994 +-006 4 8 

3-9246 - -006 4 

"8 

90° I -oooo 4-0000 2-149 I '075 '93° +'075 562 1052 

4-6190 — -070 490 

1052 
Af=o-oo4i6 by Weddle's Rule 
=0-0040 by formula ^8' sin x- 
/i = 0-069 ^s against 0-060 in the minimum error system. 
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^=maximum positive scale value near mid-parallel, 8 = co-latitude of that 
parallel of maximum scale error, then if 

8 = -|.(j8-a), 77= (sec 8)4= 1-075, i^ = cosx = l, cos 8^ = ;^^^ = 0-5773,. 
8,^ = 54° 44'; r„ = (i- cos S/A)/«Zr= 0-7868, ^0^ = 0-6190. 

The maximum scale error has been reduced from 10-3 per cent, to- 
7-5 per cent, at the expense of increasing the. mean square scale error from 
6-0 per cent, to 6-9 per cent. It will be observed however that this 
maximum error occurs at the very edges of the map and then rapidly 
diminishes, whereas in Tissot's system the errors between the standard 
parallels are much larger, being 50 °/^ greater at the mid-parallel. 

When the mid-parallel in Albers' Projection becomes the Equator the 

formula changes into r — r^ = —,siTi <^, where (f> = go° -6 or angle from 

n 

equator. The minimum value of M will be found by taking 
«' = sin 8/8 = I — — H- — 



Map of the World in three parts in Minimum Error Conical 
System as proposed by M. Collignon. 

According to Germain (p. 103) M. Collignon proposed to construct a 
map of the World in three parts by projecting the portion of the surface 
from the poles to latitude 30.° on the Isospheric Stenoteric System of 
Lambert, i.e. the Conical Equal Area System or r= 2 sin \6jjn, and the 
middle zone covering 30° on each side of the Equator on the Isocylindric 
System. The greatest error of distance is said to be 15 per cent, and of 
angle 8° 12'. 

If we use the Minimum Error Equidistant Conical System instead, we 
have for the parallel scale error at 60° from the pole 

4-/3712= (cm. of 6o°)7i2 = 9 per cent, 

and for the middle part we have in Murdoch's Third Projection parallel 
scale error at 30° from Equator = + 8^3 = (cm. of 3o°)73, i.e. the same, so 
that the sections would have a rolling fit along the boundaries. The 

greatest alteration of angle would be 2 sin~^ — ^ = about :;°. The mean 
^ ° 2-09 ^ 

square error m the conical parts =j — -^ = 0-91. 

This is probably the best map of the World that can be obtained in 
three separate sheets. The middle portion of the map covers exactly half 
the total surface of the sphere and the end portions each one quarter. The 
scale along the Equator is of course too small by half the excess at latitude 
30°, or by about 4-J per cent. 
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G. W. Hill's Projection. 

The celebrated American astronomer and niathematician, G. W. Hill, 
in a Paper published in the Annals of Mathematics, Second Series, Vol. lo. 
No. I, October iqo8, investigates a projection which he claims to be the 
best possible on Tchebycheff's (or Chebychev's) Principle that "A map is 
the best when, within its limits, the logarithm of the scale has the least 
possible oscillation." This Paper is rather difficult, but it appears to be 
stipulated that the map shall be conformal, i.e. orthomorphic, and the 
principle leads to the result that the scale error must vanish all round the 
boundary. 

The scale along the parallels is no longer uniform and the meridians in 
addition to having a variable scale are no longer straight lines. 

Hill applies his results to an example extending from latitude 40° to 
latitude 50° and for a difference of longitude extending 5° on each side of 
a central meridian. He gives a table of scale errors of which the sum is 
o"o754 for 66 values, the arithmetic mean of which is o'ooii4 and the 
maximum o'oo254. 

If Airy's Zenithal System had been applied to the same map it would 
be comprised in a radius j8 of about 6° '3, and so the mean square scale 
error /* would be ff'l Ji<)2=o'ooo()o(), the average error = 07979/A = 0-00072, 
and the maximum error /878 = o-ooi8. 

If the Conical System were applied to the same map a zone 10° wide 
would cover it, the mean square scale error would be 8767 (where 8 = cm. 
of 5°) = coon, and the maximum error would be 873 = o-oo25. 

It does not appear that Hill's Projection has any advantage over those 
we have investigated and it is certainly more difficult to compute. 

It still remains an open question, perhaps, as to whether there is a class 
of projections, i.e. of developable surfaces, better than the Conical. The 
writer's impression is that there is not, but so far he has not been able to 
prove this absolutely. 

Best projection for a given country. 

As we have already stated, the best projection for a given country 
depends on the shape of the country. Tissot has developed a method of 
finding such a projection on the principle of making the maximum scale 
error as small as possible. He applies this to finding, the best projection 
for a map of Spain and (p. 69) he gives the resulting scale errors. The 
map extends from latitude 36° o to latitude 43° 50'. The representation 
appears to be nearly conformal. The maximum scale error is +0-00119 
and the mean square scale error works out at o-ooii6. 

If the Zenithal System had been used it appears that Spain could be 
included in a radius /3 of about 5° = 0-087 cm. If the Equidistant Minimum 
Error Projection were used the maximum scale error would be ^78=o'oo°95> 
and the mean square scale error would be ^7i3'i =0-00058. This is there- 

y. 4 
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fore better than Tissot's, for we could have conformality with the same 
maximum error and a mean square error oi fi^j J<)6 = o-ooo']']. 

If however Murdoch's Third Projection were used we could take 
/8 — a = 8° and 8 = 4° = o'o7o cm. Hence the maximum scale error would 
be SY3 = o'ooi6 and the mean square scale error would be §^^45 =o'ooo73. 
This could be made conformal like Tissot's with the same maximum error 
and by increasing the mean square error to o'ooo73 x J 2 = o'ooio2. 

Either of these projections is therefore at least as good as Tissot's and 
both are easier to work out. 

Best projection for Italy. 

Mr J. I. Craig, of the Egyptian Survey Department, in his pamphlet on 
Map Projections, raises the question of the best projection for Italy, a long 
narrow country with a general bearing of about 45°. It would appear that 
a zone about 5^° wide would cover Italy proper, excluding Sicily. If we 
used an oblique conical system we could represent Italy on Murdoch's 
Third Projection with a maximum scale error of (o'o48)73 = 0-00077 ^nd a 
mean square scale error of o-oo23o/i3'i =o'oooi7. To represent the same 
country on the Equidistant Zenithal System we should have to take /3=about 
$\° and then the maximum error would be ;878 = o'ooii5 and the mean 
square error would be /S^iS'i = o'ooo7o. 

It appears therefore that we shall only do better in the Zenithal System, 
for countries of moderate extent, if y87i3'i is smaller than 876'7, that is if 
/3 : S is less than i'4 : i, or if the country to be represented can be included 
in a parallelogram whose longer to shorter side is less than about i"4 to i. 
It must be remembered that the Conical is always, theoretically speaking, 
better than the Zenithal System, but we are excluding the conical projection 
up to pole of projection, where such pole would occur in the middle of a 
country of moderate size and produce a fan-shaped map with a small gap 
= 36o°(i-^)=iS''^l 

As we have remarked Hammer has published Tables for Conical Pro- 
jections with pole in any position. The formulae are capable of a little 
simplification when the axis is inclined at 45° to the meridian as in Italy, 
but even then this skew system cannot be recommended to cartographers, 
as the decrease of scale errors will be more than counter-balanced by the 
disadvantages of complex formulae and contrary flexure of the projected 
meridians and parallels. 

Summary of Result^. 

The Minimum Error Projections of General Conical Class have been 
thoroughly examined by investigating mathematically the general function 
as found by the Calculus of Variations. 

As in the Zenithal Class it has been found possible to simplify this 
formula by expansion and by the stipulation of making the total area of 
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the map true, with no increment in the mean square error to the first or 
even second term of its expansion. 

The Minimum Error Conical Projection of a zone is then found to be 
identical with the Third Projection of Murdoch invented in 1758. 

This can be made Orthomorphic or Equal Area as may be desired, 
with no change in the depth of map or radius of mid-parallel, by small and 
easily calculated changes in the radii of the intermediate parallels, thus 
producing two standard parallels at about one-fifth the depth of the map 
from the limiting parallels. 

This class includes the Cylindrical Projections, Direct and Transverse, 
in which the property of Minimum Error will also be produced by making 
the depth of the map true. 

When the Conical Projections extend up to the pole the principle of 
making the total area true also apphes, but the Equidistant Projection is 
not quite so near the absolute Minimum Error as in the zone map. It is 
however still retommended as the best in- practice on account of its 
simplicity. 

The general practice hitherto has been to make the maximum error in 
this class as small as possible. This practice leads to a maximum error 
one-twelfth smaller than in the system now investigated, but the average 
error is 15 per cent, greater, and the valuable property of truth of radial 
scale and total area is lost. 



The best projection for a Cadastral Survey. 

A good many of the refinements we have been considering may appear 
unimportant and unnecessary to the practical map-maker, so long as his 
work must be plotted on a substance Hke paper which changes owing to 
hygrometric conditions in a way that cannot be allowed for, from day to 
day, and even from hour to hour. It is prudent however to carry our refine- 
ments to a degree further than the limits of such changes, and we should 
always look forward to the possibility of some substance less changeable 
than paper, as we know it, being invented. 

But there is one case in which the greatest care in selecting the best 
projection should be exercised and that is the projection on which to base 
the rectangular plane coordinates of a cadastral survey. 

In the great self-governing Colonies and in most of the Crown Colonies 
the system of land tenure is based on accurate survey. Boundaries of pro- 
perty are defined by the bearing and distance of lines between fiducial 
points on permanent marks set at intervals round the property. The survey 
round this boundary is closed to a prescribed limit of error by the method 
of balancing latitude and departure, and one or two points of it are con- 
nected to the general triangulation of the country either directly or by 
intermediate traverse survey to a trigonometrically fixed point. 

4—2 
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The Cadastral Survey Department does not want to be worried with 
questions of spherical excess or excentricity of the Earth, and in supplying 
it with a system of plane coordinates we ought to ensure that the errors 
of the projection chosen will lie well within the limits of accidental error 
allowed in the Cadastral Branch. 

An illustration from the Survey on which the writer was engaged may 
not be out of place. This was in a tropical country for the most part 
undeveloped and covered with virgin forest, even up to heights of 5000 feet, 
consisting of huge trees from 100 to nearly 200 feet high, intertwined near 
the ground by a dense mass of undergrowth and creepers, through which 
it was impossible to force a way unless preceded by a man with a clearing 
knife. The country was covered with a network of major triangles by the 
Trigonometrical Branch, of sides averaging say 20 miles in length, and 
very much of the interlaced type used in Great Britain and Ireland. The 
clearing of hill-tops and suitable points on the lower ground for this purpose 
involved great expenditure of money and labour. This major triangulation 
was however of very little direct use to the Cadastral (or, as it was called. 
Revenue) Branch, as settlement in such countries always takes place at 
the beginning along rivers, roads and railways, which are necessarily on 
low ground. A general system of minor triangulation was in this case ruled 
out by the great cost of clearing. But, owing to the development of the 
system of accurate traversing with the long tape of steel or invar, first in- 
stituted in Australia, the difficulty was overcome. The method used was 
to run high class traverses, locally called standard traverses, along roads, 
railways, and rivers (and sometimes watersheds), leaving permanent marks 
in pairs or threes every half-mile or so. The traverses were connected with 
the major triangulation at intervals where suitable opportunities occurred, 
generally frorn 7 to 10 miles apart, either directly or by the fixation of 
suitable minor points. The limit of error allowed in such standard traverses 
when closed on trigonometrically fixed points was i in 10,000, but the 
average yearly limit attained was near i in 30,000. These traverses then 
afforded a frame- work ready at hand in the most convenient places for the 
Cadastral Branch, with which to connect their boundary surveys. The hmit 
of error allowed in the Cadastral Branch was i in 4000, but the average 
attained was more like i in 10,000. 

The projection used in this survey was the same as that used in Great 
Britain, India, South Africa, etc., namely the Cassini or Transverse Cylin- 
drical Equidistant, cast on an average value of the Earth's radius in the 
vicinity, and with only^one standard (so-called) parallel, i.e. the meridian. 
We have seen that in this projection the x or N-S coordinates are always 
too large in the ratio sec B : i, where 6 is the spherical angle of a point 
away from the central meridian. When 6 is small the error is therefore ^6^. 
If 6 corresponds with a distance of 50 miles, as it did in one instance in 

this Survey, this error amounts to -{7-^ ) = i in 12,300 nearly, which 
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is rather uncomfortably near the limits of error in the traverse surveys. All 
areas on which land-tax is paid are increased to the same degree, which is 
an inconvenience to the lessee at any rate. Of course this error was allowed 
for in stating the closing errors of the traverse surveys, but the smaller it 
can be made the better, and our investigations have shewn that by reducing 
all N-S coordinates in the proportion ^8^, where S is the limiting spherical 
angle of the projection from the central meridian, the maximum error is 
reduced in the ratio 2 : 3, and the total area of the projection at any rate 
is made correct. 

It is curious that this device of making standard two parallels (or sub- 
meridians) in the Conical or Cylindrical Projections when used in national 
surveys, does not seem to have been noticed before. The only mention the 
writer can find of the use of such a device is by Jordan in referring to the 
Survey of Mecklenburg (p. 350) where the Conical Orthomorphic Projection 
with initially one standard parallel was used, bi^t by a little device (durch 
einen kleinen Kunstgriff), which was evidently the use of a scale constant 

I — , the maximum scale error was reduced one-half 
4 

The Cadastral Branch, as stated, was supplied with the rectangular 
coordinates of the trigonometrically fixed points and the points fixed in 
the standard traverses, in this Cassini System, but in supplying bearings 
from triangulation stations such bearings were given from a meridian at 
right angles to the great circle from such station to the principal meridian, 
according to the practice usually followed hitherto. Thus the bearing from 
a point P to 2, point would differ from the bearing from (2 to /" by the 
spherical excess of the space pPQq, where p and q are the feet of the 
perpendiculars from P and to the principal meridian. This difference 
might amount to 13 seconds of arc for a distance PQ= 20 miles running 
parallel to the principal meridian and 50 miles east or west of it. This is a 
small amount compared with the accumulated error of a traverse 20 miles 
long, but still such diiferences are rather puzzling and inconvenient to the 
Cadastral Surveyor, and having once chosen a projection it would be better, 
in the writer's opinion, to supply all bearings as well as coordinates, on the 
plane system. 

For a country nearly round or square, from which the triangulation is 
not likely to be continued in the N-S or E-W direction for any purpose to 
neighbouring countries, the Minimum Error Equidistant Zenithal Projec- 
tion has many advantages. As already stated the maximum error would be 
less than in the Cassini in the ratio 8 : 3 and the amount of spherical excess 
to be forcibly dispersed would be about one-half. The projection is easily 
calculated for it is only necessary to calculate by spherical trigonometry 
the bearings and distances of all points from the centre of projection, and 
then to reduce such distances in the ratio i - st/^^ ■ i> which means sub- 
tracting a constant from the log distances, before working the coordinates 
by simple sine and cosine of the bearing. 
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Projection for Surveys of Hanover and Egypt. 

As already remarked in the surveys of Hanover and of Egypt the 
Cassini Projection was rejected in favour of the Conformal, or Orthomorphic,, 
Cylindrical System of Gauss. This differs from the Cassini in the fact that, 
in order to produce the same scale error in the E-W direction as in the N-S 
direction, the coordinate y in the former direction instead of being the true 
distance on the sphere is increased in the ratio 

/"seced'^=6l(i+ie2+-i,V^*...) 

to Q, that is y is increased tojfn-^^...], where R is the Earth's radius 

in the vicinity. The projection as used in the surveys mentioned has only 
one standard parallel, i.e. the principal meridian. It could be improved in 
the same way as the Cassirli by multiplying both coordinates by the scale 
constant (i —\'^\ where S is the limiting spherical angle from the principal 
meridian, which would also make the total area correct. 

In the Orthomorphic Projections the errors of bearing for small distances 
round a point are practically zero and for large distances are much smaller 
and easier to compute than in the Equidistant Projections. The errors of 
distances are also easier to compute. For these reasons the writer is 
inclined to think that for Cadastral Survey purposes these advantages more 
than counterbalance the extra initial labour of computing the radial co- 
ordinates, and that Surveys about to start a new or revised system of 
projection would be well advised to consider the merits of the orthomorphic 
projections, whether direct or transverse, and to adopt the scale constant 
which renders the total area true and the average error a minimum. The 
advantages of the orthomorphic system have also become apparent during 
the war for artillery and aerial purposes, the new projection for the war 
maps of France being based on the Conical Orthomorphic System with 
two standard parallels, though with Tissot's criterion of least maximum 
error instead of the criterion of least average error as recommended in this 
Paper. This map being more liable to extension in the E-W direction is on 
the Direct System and the standard parallels are actual parallels of latitude. 

The Finite Errors in Conical and Cylindrical Projections. 

So far we have only discussed the differential scale errors at a point, 
but as Hinks has pointed out (Chapter IX) these are perhaps not of much 
use in practice, when the question mostly arises of finding the actual errors 
in distance and bearing between two points at finite distances apart on the 
map, as compared with the corresponding quantities on the sphere or 
spheroid. Formulae for such finite displacements have been investigated 
to several terms of their expansions in the works of Jordan and Helmert 
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for the Cassini Projection, and in Jordan for the Gauss Conformal (or 
Transverse Cylindrical Orthomorphic), and Direct Conical Orthomorphic 
Projections. Jordan also investigates the coordinates and finite errors of 
the Direct Cylindrical Projection, Equidistant and Orthomorphic, tangent 
along the prime vertical great circle at a point, and called by him the 
"Querachsige Projektion." This, as he shews, only differs from the Conical 
Projection, tangent along a parallel of latitude, by small quantities depend- 
ing on the fourth powers of the coordinates. The latter has one advantage 
over the former in that the projected meridians are absolutely straight lines. 
In all these investigations only one standard parallel is assumed, but the 
results will apply to the case of two standard parallels without change for 
the Orthomorphic Projections so far as the finite angular corrections are 
concerned, for the scale constant is applied to both coordinates, and so far 
as distances are concerned it will simply be necessary to deduct the scale 
constant from the distance ratio as given by the formula of Jordan. In the 
case of the Equidistant Projections, Direct and Transverse, the formula for 
two standard parallels will become more complex as will be shewn later, 
but the scale constant will tend to reduce the amounts of the corrections. 



A mathematical paradox in the General Formula of 
Minimum Error Conical Projections. 

There is a rather curious mathematical paradox in the general formula 
for the Minimum Error Conical Projections which deserves notice for its 
mathematical interest, though it is not a question that need trouble the 
practical map-maker. It will have been noticed that in the solution of the 
general case of a map covering a zone from a lower limit a to an upper 
limit /8, the Calculus of Variations exacts that the radial scale error shall 
vanish at both limits, i.e. drjM= i both when 6 = a. and 6 = ^. If however 
we make a = o we approach the case of map extending to centre and in the 
case of the Zenithal Projection, as amended by Clarke, we get for the radial 
scale at centre, not the value i, but the same as the perpendicular scale 

Wsin 6 = A = i-^.... 

' 10 

Again, in the General Conical Class the condition of Minimum Error 

leads to two standard parallels and a mean square error = —j— ; but if we 

start with the assumption that the map extends to the centre, we are led to 
a projection with one standard parallel only and a mean square error 



_ . Now if in the former value for \i. we allow a to approach the 

\/768 

value zero, then S approaches the value \^ and —, — approaches the value 

—^ . which is slightly larger. Professor Nicholson was kind enough to 

V720 
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give some attention to this paradox, and without going into a full investi- 
gation, the explanation seems to be that so long as we have a value fer a, 
however small, the projection will have the radial scale true at a and two 
standard parallels ; but when a is actually zero a discontinuity takes place 
and we have a projection with only one standard parallel. The paradox 
seems to be somewhat similar to that of the deflection of a circular elastic 
plate supported round the edge. If there is a hole in the centre however 
small we get a different result from that derived from the initial assumption 
of plate continuous to the centre. 



CHAPTER III 

THE CONICAL ORTHOMORPHIC PROJECTION WITH TWO 
STANDARD PARALLELS (LAMBERTS SECOND) FOR THE 
SPHEROID 

We have hitherto considered the Earth to be a sphere of radius R, taken 
as unity in our formulae for brevity, but as it is usual, especially in the 
projections employed for the rectangular coordinates of triangulation, to 
take account of the excentricity, we will now investigate the expansions of 
the spheroidal formulae for the above. The expansions of the formulae for 
the case of one standard parallel are given by Jordan (pp. 521-4), Pizzetti 
(pp. 389-gi) and Oscar S. Adams, Special Pubhcation No. 47 U. S. Coast 
and Geodetic Survey. 

The formula for the spheroid in terms of our original notation is 



'-^<-*»K^^^f. 



where e is the Earth's excentricity. It will be now more convenient to put 
.^for the constant 2 A and write 



^^<'»*»K^^fSir 



Calling the radius of curvature to the meridian p, and the grand normal or 
radius of curvature at right angles to the meridian v, we have ^ p' 

p = «(i-«^)/(i-^cos=e)*, and v = a/(i-«2cos^6l)4, 

where a is the semi-axis major of the spheroid. The scale along the meridian 
or what we have hitherto called the radial scale, is given by drjdm, where 
dm is an element of the projected meridian. Now 

J ij JIM jaij ^ ^^ dr {1 - e^ cos" ef 

dr dm = dr do . dd dm = - -^ = — - ^ '- . 

' ' ' p do dd a(x—e^) 

Differentiating logarithmically we have,putting r'iox drjdQ, r"{or d^r/dd'', etc., 

fl€ 

\ogr = \ogK+ ;^ log tan ^^ + — {log(i + «cos 6) -log(i-fcos6)}. 



r /I r sm \ _ ra ( ; 

r ~ \sin 6 i -e^ cos^ 6/ sin 6 (i - e^ cos^ 0) ' 

(i-e^cos^6i)^ 



Therefore drjdm = nr - 



a sin 6 



and as the projection is orthomorphic the scale along the parallels has the 
same value. Putting as before p( for the co-latitude of the mid-parallel and 
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S for the angle from x to standard parallels ^i, 6^, we have 0-i_=x-K 
02 = x + 8. The condition that the scale at 6^ shall be the same as the scale 
at O2 is 

(i-e^cos^^ii)* (i-e^cos^(92)4 

a sm 61 a sm 0^ 

where r-^, r^ are the radii at 61, B^. Writing out the full expressions for 
i^i, ^2 and cancelling common factors this leads to 

log sin (9^ - log sin ^i^ - log ( i - e" cos' Qj)^ + log ( i - g^ cos^ &^^ 

n J ■ 

,. , 1^ , /I +« cos 6*2X2 . /i + gcos6lA2 

logtani^2-logtanie, + log(-^^-^^^J-log(^^^^^J 

The additional condition that the scale shall have the value unity along 
these parallels is 

( 1 - «= cos^ ^1)* ( I - fi' cos' 6^^ 

a sm 61 a sm 0^ 

which leads to 

ne lie 

„ Visin^i . fi—ecos6i\i Cgsinl 

li. = - 






n (tan J6'i)'' \i + « cos Oj n (tan ^^2)"' \i + ^ cos ( 

It will also be noted that the radii of the standard parallels have the simple 
values ^1 = Vi sin Q-^jn, r^ = v^ sin 0^\n. Such are the complete formulae for 
two standard parallels chosen beforehand. To calculate the values of r by 
them with exactitude log tables of at least 8 figures are required, so that in 
practice it is more convenient and even more accurate to use expansions, 
involving the differences and values at the mid-parallel. These expansions 
can be made in terms of either the difference of latitude A^ following 
Jordan or the difference of meridian distance Aw following Pizzetti, the 
U.S. Survey and others. These writers consider the case of one standard 
parallel only, where n = cos \, 

ne 

j^^ Vx""X / i-gcosx \ 2 
« (tan Ix)" Vi + « cos x/ 

It may "be rernarked that if we write 

. ,y . ,fl/l+^COS 

tan|4 = tan J» 



\I — (?COS ( 

then ^ is very nearly the geocentric co-latitude 

= e -f Jg' sin 26I -1- 1£^ sin 26 cos' 61 + . . ., 

whereas by expanding the above expression we find 

i=e+ |e'sin26i-H,;\«^sin2ecos'6i+ ... 

as pointed out by Herz (p. 246). If tables of geocentric co-latitude are 
used therefore which take account of the term in e* the error is very small, 
the maximum being only o"'5 for 6 = 30°. 
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By expanding the expressions involving e we can write 

_ logiosingg-logipsin 6-^^-\M^{co%^Qi- cos^^i) + \Me'^(c.Q^^O.i- cos*gi)+ ■ . ■ 
logiotanl^a-logiotan J^i+ ^e- (cos ^2- c^os ^1) + ^7l/«* (cos^^a— cos'' ^1) + . . . ' 

which is fairly easy to compute, M being the modulus of common log- 
arithms. 

But for our purpose we require to expand the whole expression for n in 
terms of S^ and e^ which can be done as follows : 

log sin ^2 — log sin Q-^ = / cot 6d$ 

= 28 cot X 1+-^-^-+ ^j-^+... 

log tan J ^2 — log tan J ^1 = I cosec OdO 

. / S^i + cos^'v 8* ■; + iScos^Y + cos^Y \ 
= 2Scosecx 1+^ — ^^-^ + ^ -+ •••) 

rh f 28^ \ 

\^{Q.o^Q^-co^^^=-r \ sin^cos^^i9=-2Se^( sinxcosx sinxcosx ••• ) 

e^ (cos ^2 - cos 61) =— (?^ I sin 6d6 = — 28fM sin x - 7^ sin x • • • j ■ 

We shall not carry the expansions beyond 8* and e^^, but if desired 
the further proceedings are obvious. Dividing we find 

/ 8^ I 8* 2 + cos= X 2 • o 2 „ \ 
•^V 3Sin'x 15 sin^x 3 / 

"' ;; 7 8^TTcos% 8* 5 + iScos^x + eos-'x 

28 cosec X I + 7 — ^~5 1 ^~i 

'^ V 6 sm^x ''^o sin'x 

— e^ sin^ X + r ^^'^ sin^ X • • • ) 

f 8% 2 • o x^ S* i6 + 7sin^x 1 

= cos Y -^i + t(i + 4« sin-x) + — ^- ^-s •••f • 

-^ ( 6 ^ 360 sm-x J 

We shall now continue the expansion of r, but it will be convenient to 

do this in terms of Aot, the meridian distance, instead of A0, the difference 

of co-latitude. We have found already drldm = nr {1 -e" cos^O)^ ja sin 6, 

and differentiating further logarithmically, dropping the factor a pro tem., 

we have 

/' /r ^ _ -ncos6{i -e^) 

~^~\rJ ~ sin^6l(i-e^cos^^)* ' 
and multiplying by dQldm, 

fr" ir'\}\ -n cos 6{\-e^ cos^ &) 



(7 



sin^^ 



therefore (-^) 



r"\ _n{n~ cos B) (i - e^ ZQ%^ 6) 



sin^e 
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Continuing 

r" _rW__ ji-e^'cos'e (i-<?')2cos6l(^-cos6l) ) 
V~'^~"\ sO sin^6* J' 

/r^ _ r^\ _ ( (i-e^cos^6)i _ 2 cos 6 {n - cos 6) (i - e" cos^ 6)^ \ 

therefore 

('^—\ - f(i -g'cos ^^)* (« - cos g) (« - 2 cos 6){i-e' cos^ 61)* ] 

Wm~''i'siir^(7^«2) "^ sin^e r 

Again 

^iv r"r' _ f 5.;^ ( I - g' cos' g)^ cos Q sin' 61 - cos g (i - g' cos' B)^ 

T'^^^^X sin'e(i-e') 

{n— 2 cos ^) ( I — e' cos' 6)^ 2 (n- cos 

3e' cos ^ (i — e' cos' 6)2 (« — cos ff) {n— 2 cos ^) 

3 cos ^ (« — COS ff) (n — 2 cos 6) (i — e' cos' 6)* 
sin* 61 
and 

/^ _ ^LZ^A - f - cos e(i-e' cos' 61)* + 5^' cos 61 sin' 61 (i - g' cos' (9)° 
W '^ /m~'^V sin'e(i-g')' 

(« — 2 cos ff){\- e' cos' 6)'' 2 (« - cos 6) (i — e' cos' 6)' 
"^ sin'e(i-g') "^ sin'6i(i-«') 

3^' cos 6 (i — e' cos' 6)' (« — cos 6) (n — 2 cos 6) 
^ sin'6»(i-«') 

3 cos (;z -COS 6){n — 2 cos S) (i — g'cos'fl)' 
sin*6i(i-g') 



(;^-2cos6')(i-g'cos'6i)^ 2(;7-cos6i)(i-g'cos'g)^ 
sin' ^ sin' 6 



}■ 



Now 



}• 

(r"'r'\ _ i({}^/_cofey (n - cos 0) (n - 2 cos g) (i - e' cos' 6)' ) 
V ^/m ~ '^ I sin' e (i ->) "^ sin" j ' 



therefore 



riv\ _ (• - cos 61 (i - g' cos' 61/ + 5^' cos 6 sin' 6^ (i - g' cos' 6) ^ 
TU " " I sin'^(i-g')' 

4 (« - cos 0) (i — g' cos' 6y 
^ sin'^(i-g') 

3g' cos 6 (n — cos 0) {n — 2 cos 6) ( i — g' cos' 5)^ 

■^ * sin'e(i-g') ■ 

3 cos (n — cos 6){n — 2 cos 6) ( i — e' cos' ^)^ 

sin* 6 (I -g') 
n(n— cos 6) (« — 2 cos 6) ( i - g' cos' ^)'l 



i 
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- cos 6 ( I - g' cos' eY(\-e'-\e^ sin' Q) 

4 (« — cos &) (i - e^ cos' &f 
^ sin'e(i -«') 

3 cos ^ (« — cos 6)(n~ 2 cos 6) ( i - g' cos' 6)' 
sm^l 



! (« — cos 6) {n — 2 cos 0) (i — «' cos' 5)' 
sin* (9 



Now if we are expanding only as far as m!^ we can put n = cos x in the 
expressions for {/" jr)^ and {r^^/r)^^, so that only the first terms within the 
brackets will be retained. Since also 

p = a{i- «')/( I - ^ cos' 61)^ v = aj{ i - e" cos' 6)*, 
we have 

pv = a' ( I - «')/( I - g' cos' 6)', pv' = a' ( I - «')/( I - «' cos' 6f. 

Thus, replacing ^ by the value x for mid-parallel, we have 



• = ^(tanix)''( ^''''°'^ V] '"''^'~'^'°'^^^^ -' 
\ 2'^/ \i-ecosx/ t asmx 



n(n — cos x){i—e' cos' x) »2^ n{\-^ cos' x)^ '«^ 
"*" a' sin' X "2 a^sinx(i -«^) ' 6 

w cos X ( I - g^ cos' x)^" ( I - g' - 4^" sin' x) ■ni_ 



'24"'7 ' 



a''sin'x(i -^T 

We shall subsequently expand ^ in terms of x and S but for the moment 
we put 

3S6.\2 



V2 sm 6/2 / 1 ~ ^ cos c^2\ 
~ ^7tan"|^ \i +ecos6'J 

^ ^ inx(cot|xr /ijilco^XN ^ (, ^ ^,,^3 ;„ g). 
«(i-g'cos'x)^Vi+^cosx/ 

We then have radius at mid-parallel 

r, = e.'^lX (j + terms in 8), 

" «(i-e'cos'x)* 
and 

(r-cosx)(i-g'cos'x)* ^ (i-f'cos' x)^ ^' 
asinx '2 a'(i-/) ' 6 

c otx(i-^^cos'x)- / _ 4^!HJX) ^\..} (I + terms in 8) 



(;z - cos x) ^ , ^ 



= ^»? + -^^ ; ■ -- + 






V sm X 2 pv 

C( 

pv 



-^(■-4^)S-}<-^'""^'"*'- 
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Por n = cos x this reduces to the expression found by Pizzetti for one 

standard parallel, when terms in 3 vanish. Siubstituting however the 

approximate value of 

f 8^ • 1 

n = cos xii+^(i +4is^ sin^ x) • • • f 

we find the second term gives a term in w^S^ and 

(■ I w' cotx / o • 2 .m^^ 

■* (. pv 6 V ^ '^' \2 

_52i2^(i_4£2sin2;^)— ...)(i + terms in S). 

We must now expand the expression for K in order to find the terms 
in 8. We have 

fie 

g sin 61a (cot 1 6>)'' / \-ecQ%d^ <7^ 
~ « (i - fi' cos' dlj)* \i + « cos t'J 

I + cot X • 8 +7- cot X • • • ) > 



Now 
also 



;?^+ 3;>z^cosx +^2 + ^cos^X ^' 



)■• 



sin^ X ' 6 

and putting in the proper value of n this becomes 

(cot 1612)" = (cot Ix)" (i -cotx- 8 + cot^x- S' 

/2 +4C0S^X , 2 • 2 \S^ \ 

Multiplying this by the expansion for sin B^ we find 

/ 8^ \ 

sin ^2 (cot \Q^-^ = sin x (cot ^x)" ( i f «" sin x cos . 8' . . . j . 

Next take {1 —e^ cos' ^2) "^ =/say, we have 

1 r 11 / 2 2a\ ^'/^ «'sin^cos6 2 ■ a a 

log/= - 1 log ( I - «' cos' 6^), f lf= - ^_^^,Q32g = - « sm 61 cos (9, 

neglecting terms above e' ; "T ~ ( T' ) = ~ ^^ (cos' 6 - sin' ^) =-^-^ to same 
order, and -77 = 4^' sin 6 cos 6 also to same order, therefore 

(i -«'cos'^2)~^ = (i — «'cos'x)~^ (i -«'sinxcosx-8 

^/ 2 ■ 2 x^' 2 • 28^ \ 

— r (cos'^x-sm'^x)- + «'sinx cos x- -7- ••■ )• 
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In the same way we shall find 



/i-ecos^2\2 /i-«cosv\2/ „ . „ 82 

Vi+^cos^J = ( i+.cosx ) V "" ""' ''" ^^ + "'' '^^^ ^ • 2 

- ;»2e-* sm X . 7" ... I 

ne 

(I — ^ COS y\ ^ / S^ \ 

i+gcosy j ( ' + ^' ^"^ X cos X ■ S + e^ cos' X • - + terms in 8^ upwards) . 

Multiplying out we get 

ne 

/I-^COSX\% 2 2 X-i/ 2 • 2 S' 2 • 28' \ 

"iiTTcoTx/ ('~^''°^^) 4(^i + fi'sm'x- j + «'smxcosx.y + ...J, 
and multiplying this by the expression above we get 

ne 

^^ ^sinx(cotixr / i-^cosx y^^_,g,(^_^^i^ 

«(i-e'cos'x)* ^^+*'^°^^^ . 

+ terms in 8* upwards}. 

We can test the accuracy of these expansions by seeing if the scale 
works out correct on the standard parallels ; we have 

d , , f I »i' cot X / 2 • 2 \ »28^ 

dm^ '^' '^ pv 2 V '* 6 

-^^(.-4^sin^x)^'...}{-J(-^sin'x)...} 

f {i-e^ cos'' xf m^ (i—^ cos^ x)^ ( 1 + 4iS^ sin' x) cot x fnS^ 
~ \ a^ (i —e^) ' 2 a '6 

{ i-e' cos^ x)* (i - 4^ sin' x) cot x m' \ L - - Ct -^^.in^ v1 I 
d'{i-e') •6---JI 2^ «sinx;...|- 

Now »? = pS (i-- -^'sinxcosx- 8 ... j 

i2(i-e')S , . „ 

,(i -f«-smxcosx-o ...); 



(i — e^cos'x)^ 



m' = , ^ ^ f-^3(i-3«'sinxcosx-8 ...). 

(i-rcos'x) 

Substituting these and neglecting powers of e above e^ we get for the scale 
|i + (i - ^ sin' x) I' - ^ sin X cos X ■ 8' (I - i - i) . . . } 

X |i (i-f'sin'x)-"[ = i. 

which proves our formulae. 
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If we write m^ for the mean of the meridional distances from mid-parallel 
to standard parallels we have wz^ = pS + terms in e'^ upwards, so that 

8^ = — 5- to same order and 
P 

P (, - ^ .in^ v^ - <(i-g'si»x)(i-^cos'x)'' _ < 
^ «smx;- a'ii-ej pv 

up to terms in e^ inclusive, so that we can write 

(r-r,)=|.« + --^^,cotx(i-4^^sm^x)+ ^^(i + 4^^sm x).-} 

xjx-^...l. 

I 2pv ) 



The value of r-^ is 



a smx 



= i'tanx ji -y (n-ie'sin^x)}■ 



«(I — e^ cos^x)'^ 

3 

The French in using these expansions for their new maps employ up 
to the term m^jdpv only, and for this purpose it is sufficient to take for n 
the value cosx as for bne standard parallel, the correction to make two 
parallels standard being effected by the multiplier (i - m^l2pv). For r^ 
they use the value ;'tanx(i —m^J2pv). By some oversight they employ v^ 
instead of pv in their formulae. These small errors by a fortunate con- 
currence of opposing effects have an almost imperceptible result on their 
tables of coordinates. 

The writer has expanded the formulae still further as for the sphere, 
that is, neglecting the terms depending on e. For r-r„ there should be 

, , , m^ K— cos X in^ 8^ cos^ y , , , . ^ 

added -1- ^ — r-5 —52 ■ „ , and the reducmg factor becomes 

i20i?^ sm^X 36i?^sm^X 
§2 . S'' I - 2 cos^ X 

I (i - e^'sin^x) + ^^ ■■■'> the further term in n has already 

2 ^ '^^ 24 sm^x 

been given. If we take x = 9°° ^^^ V^^ sinx=i, cosx=o we obtain the 
formulae for the Mercator Projection for the spheroid with two standard 
parallels in latitudes + S. 



Minimum Error formulae for the Spheroid. 

In the case of the minimum error formulae for the spheroid we do not 
select the two standard parallels beforehand, but we make the scale values 
equal on the upper and lower parallels for one condition determining n, 
and for the other determining K we make the total depth of the map true. 
Now we have already in the above investigation determined the formula 
for n, because we have only to consider 8 as being no longer the angle 
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from mid-parallel to standard parallels but as the angle from middle to 
■extreme parallels. To find K we have 

We have already indicated the method for expanding these expressions, 
which involve long and complicated differentiations and it will suffice here 
to give the result, 

ne 

^^ V sin X (cot Ix)" / i-gcosx y L _V!:l_ \ 
n \i+ecosx/ 1 6pv"J' 

(m ^\ 
I - ^ — ^ j , where m^ is the mean of the 

meridional distances of upper and lower parallels from the mid-parallel. 
The corresponding value of r-^ is 

~T^~J^h~^V --•(i-«^sin^x)---|=''^^"X ji -r(i+|sin'x)--} . 
«(i-(rcos^x) *■ " > 13 I 

-where 8 is now the angle from middle to, extreme parallels. 
It is obvious that by writing 



r~r^ = m — ^ and K= ...\\ + — ^ > , 



6pv 



we shall obtain a projection which is (approximately) Equal Area such as 
would result from expanding Albers' minimum error formulae for the 
.spheroid. 



CHAPTER IV 

THE POLYCONIC PROJECTIONS 

The Simple Polyconic Projection, which the French cartographers call 
the "projection polyconique ordinaire des Americains" because it seems 
to have been first introduced by the U.S. Coast and Geodetic Survey, has 
the following characteristics. A standard meridian of longitude is chosen 
which is usually the centre line of the map. This is developed in its true 
length on the map and each parallel of latitude is developed on its own 
tangent cone opened out on the plane of the map. Thus the standard 
meridian is true to scale on the map, and all the parallels cut it at right 
angles and preserve their true scale ; but because the radii of the projected 
parallels vary as the cotangent of the latitude from point to point, the 
developed meridians on each side of the standard are no longer straight 
but curved lines concave to the central meridian and their lengths are too 
great in proportion as their difference of longitude from the central meridian 
increases. When this projection is used, however, for only a small distance 
on each side of the central meridian its errors are correspondingly small, 
and it has the advantage that concise Tables can be calculated once for all 
for plotting the coordinates of the graticule of any map. Extensive Tables 
of this kind have been published by the U.S. Coast and Geodetic Survey, 
for each minute of latitude from the Equator to the Pole, and for intervals 
of longitude up to 30° from the central meridian, in the Report of that 
Survey for 1884. Full particulars of the theory as applied to the spheroid 
will be found in Special Publication No. 57 of that Survey. These Tables 
were not used by the Committee for the International Map, as different 
values of the Earth's spheroidal elements were thought desirable. It seems 
to the writer very desirable that some international agreement should be 
come to in matters of this kind, as many valuable Tables have been 
calculated in different countries using different values of these elements, 
thus seriously vitiating their usefulness for purposes of comparison. As the 
actual geoid differs considerably in places from any of the spheroids put 
forward from time to time, but as these spheroids differ among themselves 
by quantities so small as hardly to concern the map-maker, the best course 
would appear to be to decide on which of them the most work has already 
been done in the way of calculating Tables and to adopt that spheroid for 
all international purposes. 

The equivalence of the Simple Polyconic Projection with the Cassini 
will be readily seen by comparing the expressions for the expansions of the 
coordinates. Those for the Cassini are given by Jordan (p. 483). The 
origin is supposed to be at the Equator, x is the coordinate in the NS and 
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y the coordinate in the E-W direction positive N and E, <^ is the latitude 
and X the longitude (reckoned from the central meridian in cm.) of a point 
on the map, M is the meridian distance of the point from the Equator 
(called B by Jordan), v the grand normal at latitude ^. Then we have 

Cassini x = M+ | vA.^ sin ^ cos <^ + Jj v\* sin <^ cos <^ (5 - 6 sin^ <!>)+■■■, 

Polyconic x = M+ ^vA.^ sin <^ cos </> - ^\vA.^ sin" <^ cos <^ . . ., 

Difference -^^^vX* sin <^ cos' <^ 

Cassini jc = vX cos <^ - JvX"sin^<^ cos <f> - y^vX' sin^ <^ cos <^ (8 - 9 sin^<^). . ., 
Polyconic y = vX cos ^ — ^vX'sin^ (^ cos kJ) + yJ^vX'* sin"* cos ^ . . ., 

Difference tV"'^" ^'"^^ 4' cos' <^ 

In the International Map of the World on the scale of i /million the 
device was adopted of reducing the scale along the central meridian by 
half the amount of the error at the extreme meridians, that is, by 
^X^cos^^OT, where (j>m is the mean latitude of a sheet. This appears to 
have been regarded as a rather empirical correction by its proposers, but 
the investigations we have made prove that it was really based on a logical 
criterion, that of making the maximum error as small as possible. In the 
writer's opinion, however, it would have been better to take as the reducing 
factor JX^cos^^m, following Airy's criterion of making the average error a 
minimum and the total area true. But as a matter of fact the correction is 
on these sheets so small as to be within the errors caused by hygrometric 
changes in the paper. It may be noted here that in these Instructions and 
in Germain, x is used for the E-W and y for the N— S coordinates. The 
writer prefers to use these symbols the other way about, for if a is the 
bearing of a line s nearly always reckoned from the North in survey work, 
we have the common forms x = s cos a, y = s sin a. 

Equal Area and (approx.) Orthomorphic Polyconic Projections. 

J. I. Craig has given on p. 61 of his Theory of Map Projections the 
formulae for a Polyconic Projection which is Equal Area. Instead of the 
inclination of the projected meridians being X sin <^ as in the Simple 
Polyconic, this quantity becomes a series of which the first two terms are 
, . / vX^cos^rf> \ 

By changing the sign of the second term and taking this quantity as 
X sin ( I + ^- — 7 — ~ ■■■] we obtain a polyconic projection which is ap- 
proximately orthomorphic. It will be found to agree in fact with the 
expansion of the Transverse Mercator or Gauss Conformal Projection to 
the first term of the scale errors. 

It would appear to be impossible to obtain a truly Orthomorphic 
Polyconic Projection without sacrificing uniformity of scale along the 
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central meridian. Such a projection was invented by Lagrange and is 
usually known by his name. That it really belongs to the Polyconic Class 
is pointed out by Tissot and Herz, and also by Oscar S. Adams in Special 
Publication No. 57 of the U.S. Coast and Geodetic Survey. 

Lagrange's Projection for the Spheroid. 

The full investigation of Lagrange's Projection for the spheroid will be 
found in Fiorini (pp. 533 et seq.) and U.S. Publication No. 57, and it is 
not necessary to repeat it here. The principal meridian coordinate x can 

be written in the form x = ^ jij, , where A, k and n are constants 

Tie 

^ ) , u being 

the co-latitude of a point on the meridian. The scale value dxjdm will be 
given by znAjv sin u {kc^ +2 + t^jk), and Lagrange chose the constants k 
and n so that the first and second differential coefficients of this quantity 
should vanish at the point of co-latitude u,, chosen as the centre of the map. 

The first condition leads to ku'' = " . From the second condition 

n — cos Mq 

Lagrange derives n = (i +%\v?u^^, but in doing so he neglects the excen- 
tricity. Captain G. T. McCaw and the writer almost simultaneously found 

the correct value for n should be « = {i -1- sin^ u,, + 2e^ sin^ ^'0/(1 - ^)} • The 

constant A is usually chosen so that the scale shall be true at co-latitude u^, 

but our investigations shew that a better result would be achieved by taking 

A so that the scale error along the central meridian should be one-third 

that of opposite sign along the extreme meridians. 

For a map of small extent in longitude Lagrange's formulae, when 

expanded, agree to the first term of scale errors with those of the Transverse 

Mercator, so there is no advantage in using it for maps of moderate size. 

The variation of scale along the central meridian is very slow. The writer 

finds 

_ . / Aa' 3 — 2 sin^ «o ^u^ \ 

Ax = R!^u ( I -t- cot ?/o — ^^r^ — ' ■ . .. I , 

\ 24 sm' z^o 120 / 

so that if the centre of the map is on the Equator even the third differential 

coefficient of the scale value vanishes. 

This projection has some interesting geometrical properties, as the 

projected meridians and parallels are true circles even on the spheroid. It 

has this advantage over the Mercator that a map of the globe can be 

constructed on it without the scale error becoming infinite (which applies 

also to all the Polyconic Projections) but the exaggeration from 80° latitude 

to the Pole is enormous (compare Projections Nos. v and viii in Germain, 

Plates VI and VII). 



CHAPTER V 

THE FINITE ERRORS OF PROJECTIONS 

As we have already remarked in Chap. II, the differential scale errors at 
a point of a projection are not perhaps of much value to the users of maps. 
Generally speaking, the projection should be so chosen that the greatest 
error in distance and bearing between two points at a finite distance apart 
taken from the map shall be smaller than the accidental errors due to 
stretching of the paper, etc., and smaller than the degree of accuracy 
required in the purpose for which the measurement is required. But when 
the projection is required for the purpose of fixing points with great accuracy 
by rectangular coordinates as for a Cadastral Survey, or even for artillery 
work, it may be necessary to know how the distance and bearing derived 
from such rectangular coordinates differ from their true values on the 
surface of the sphere or spheroid. 

For the Gauss Conformal or Transverse Mercator Projection, as used 
in the Surveys of Hanover and Egypt, such formulae have been worked out 
for the spheroid to a high degree of accuracy by General Oscar Schreiber 
in his Theorie der Projektioiisniethode der Hannoverschen Landesvermessung 
(Hannover, Hahn'sche Hofbuchhandlung, 1866). They are also given as 
far as the fourth powers of coordinates by Jordan (pp. 507-12). These 
formulae apply to the case of scale true along the meridian (or what we have 
called one standard parallel in our investigation of the Conical Projections), 
but in the case of orthomorphism they can be applied to the similar projection 
with two standard parallels or sub-meridians without any change in com- 
puting the angular corrections, and for the distance correction it will only 
be necessary, after computing the correction by these formulae, to apply 
with opposite sign the scale reduction factor. 

The same formulae can be used, with simple interchange of the x co- 
ordinates for the jv, at any rate to the first term of the errors — which ought 
to be near enough for all practical purposes — to the Direct Mercator, the 
Direct Conical or Lambert's Second, and the Cylindrical Tangent to a Prime 
Vertical (Jordan's Querachsige) Projections. 

As to the Cassini, standard along one meridian, Jordan besides others 
gives the formulae for the finite errors for the sphere but not for the 
spheroid. These formulae, to the first term at any rate, can be applied to 
the spheroid by substituting pv for R^, where p and v are taken for the mean 
latitude of the two points in question. The correct formulae for the 
spheroid, as far as the fourth powers of the coordinates, are given by 
Helmert (equations (5) and (6), p. 417), after a rather difficult investigation 

S-3 
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and in a rather difficult form. The writer succeeded in checking Helmert's 
result in two ways, first by a transformation from the coordinates of the 



B(.i;,yB) 



q geodesic 




Fig. 2. 

Gauss Conformal, and secondly by pursuing the method of Gauss (see 
Schreiber, p. 39) of forming the differential equation of the curve assumed 
by the projected geodesic between the two points, by the Calculus of 
Variations. This differential equation is more complex than that in the 
Conformal System, but a solution can be found in series. It can be 
written thus: 

d^f _ y ( ^ 4{x — Xi) rft\^ 



(cos^ /i + 2 sin^ j^i cos t-^ ■ 



2rpty^ sin /i cos^ t-^ 



the symbols having the same meaning as in Jordan, p. 511, where an 
abridgement of Schreiber's method is given. This equation when solved 
will give the angles between the chord joining the two points on the 
projection and the tangents to the projected geodesic ; but to these must in 
addition be applied the' change in the bearing of this chord due to the 

1 V ^ 

projection not being orthomorphic = - ^^^ ^arAy. Finally we shall arrive 

2 7^-^ S 

at the following equations, where 

Ci.x = X2- Xi, Ay =y^ — ji'i , s' = Ax^ + Ay'*, 
5^ = square of geodesic distance between the two points on the spheroid, 
r-^ = pv at the point x^ y-^, 

^m = P" at the mean latitude of the two points x^y^, x, r,,, 
<j> = that mean latitude, 
f = tan 4>, if = fi^ cos^ <^/( I — e^). 



S,= 



A^ A .y A 1; 

(l)'i + Ajc) + ;r;7yf- (jCi' +JF1JC2 +>-/) 



6r 



C*^m ^9. 



Tft^x^y-^ fftt^x^l^f 
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A„ = s, + 8, = ^^iiZLii:?) _ '?^^^^!Ay 

= spheroidal excess of space A^ABBf^, 
^^ cos^ /i = Ax- ^ (SJfi + 6^1 Aj + 2 A_y) + 



^=sin=A = A/.^,(3;',A^.A/).^-^!^4^..., 

•S'' = Ax= + Ay - -—5 (371' + 37, Ajv + A/) + -^--^ (2^1 A7 + Ay) 
3'm 3' 

Ax^ Tj^/Ax^AjK , , 

= Ax^ + Ay- -—2 (71' +j'ij>'2 +;'2-) + , , - (jKi +y^), 

6'm 3' 

which agree with Helmert's formulae remembering that he reckons x and j' 
as positive S and W. 

The above expressions, as well as those of Helmert, Jordan and 
Schreiber, refer to the geodesic lines on the spheroid and their directions. 
Before using them the appropriate corrections should be made to the 
observed directions, or we may add these corrections to the above expressions, 
thus giving the observed directions. This amounts to adding ^rfA.xAy/r^" 
and ^rfAxAy/rJ' to the expressions for Sj and 83, i.e. instead of writing Ay 
in the first term of each we write Ajc (i +Tf). As these corrections are so 
small as to be generally neglected we may well reject the higher terms in rf, 
and use only the terms involving the second powers of the coordinates, with 
the proper value of r^. 

In the case of the Minimum Error Cassini Projection, our Ax is reduced 

by ^y/r^^Ax, where F is the limiting value of the y coordinates. The 

effect of this on the bearing of the direct line from x^y^ to x^y^ is to deduct 

AxAi'K^ 

~~— . We shall therefore have 

. Ax{2y-^+y^ Ax Ay 

. Ax(jFi + 2y,) Ax Ay 

h = g' 2 -^ ' - 67^ W +y^y^ +yi-Y')-.... 

Similarly we shall have 

^= Ax^ + Ay -, {m' + iyi^y + V- y") + ••■= 

or S = s[i~^Ay^'+yry, + y.'-Y') + ..)^, 

tending to reduce the errors in both cases. 



>2 FINITE ERRORS 

If, instead of making the average scale error a minimum, we make the 
maximum scale error at the extreme meridian Y equal and opposite to the 
scale error along the principal meridian, as in the International Map 
on the scale of i-Million, we have only to replace the V^ in the above 
formulae by fF^ to obtain the corresponding finite errors of bearing 
and distance. 

The expressions can also be written in the form 

^i = fi4l{^«(2J>'i+J'2) + (jl'i'+_ViJ'2+J)'2'- F^)sinAcosA ...}, 
^2 =T—A^x {y^ + ays) - {y^ +y^y^ +y^ - V) sin A cos /i . . . }, 

For the orthomorphic projections they reduce to 

^1 ^ 2 — 2 (2^1 + J2) + terms in 4th powers of coordinates, 



-^ = ^ {i - 6^ W +y^y^ +y' -y')-]. 



remembering that Sj and Sj will be in circular measure, and that if we want 
the scale error along the principal meridian equal and opposite to that along 
the extreme meridian we must write f V^ instead of Y^. 

Pizzetti on p. 391 gives formulae for the error in distance and bearing 
of a hne, on the Direct Conical Orthomorphic Projection, starting from the 
mid-parallel with azimuth a^, which appear to be erroneous. For 5-5 he 

f COS" Q. f^ OOS fl 

eives r-r-" which should read -— — r— ^ to agree with the above (putting 

2paN„ 6p„iV„ 

\ T^ c , ■ ■f^ cos'^ a„ 

yi = o and y^ = x^ m direct system = s cos t^). For 61 he gives e = — — r^ 

s^ sin a„ cos a. 



which should read 



6p„iV„ 



An error in Jordan's Handbuch. 

In checking the formulae above the writer had occasion to use the 
formulae for Reversion of Series given by Jordan on p. 191, and in doing 
so found an error in the third line of terms for the value of x printed thus : 
/ZBU" 6BbC 2Bc i^CD SBC 2bD ^BE > \ 

■2 CD 

The coefficient of the fourth term is in error. It should be .^ ^ instead 

°^A^- 



CHAPTER VI 

THE CONVERGENCY OF MERIDIANS 

The question of Convergency of Meridians has always presented some 
difficulty to the writer, and probably to others also, especially when applied 
to the spheroid. Jordan gives two or three definitions of it in different parts 
of Vol. Ill of his Handbuch, and yet another definition at p. 501 derived 
from Gauss when dealing with the Conformal Projection. This latter he 
distinguishes by calling it the Convergency of Map Meridians, being the 
angle between the tangent at a point of a projected meridian of the spheroid 
and the principal meridian of the map. This of course is the convergency 
which a user of the map or projection wants to know when converting 
a bearing calculated from the coordinates into a true bearing. It will be 
noted however that the convergency so derived differs from the convergency 
as derived from the Cassini Projection by a small quantity =fX^»?^ sin <^ cos^ <^. 
An explanation of this difference is given by Jordan on pp. 503-4 which 
for some time puzzled the writer and it was not until he had studied the 
work of J. I. Craig on the theory of the indicatrix that the explanation 
became clear. 

It might appear that the x coordinate in the Gauss Conformal Projection 
was the same as in the Cassini. It is so in the case of the sphere, but not 
in the case of the spheroid. By expanding the formulae for the coordinates 
on the Cassini after the manner of Jordan still further the writer finds 

x = M+^vX^sm<pcos(t) + ^^v\*sin(l)COs^(t>{5 -^^ + 5v')+ ■■■> 
y = vX cos <^ — ^v\^ sin^ <^ cos <^ — ti^-^vX^ sin^ <^ cos* ^ (8 - ^^ + Srf) — ..., 
whereas on the Gauss Conformal we have (Jordan, p. 497) 

x = M+ ^v\^ sin <l> cos <ti + -^vX* sin <^ cos' <^ (5 - /^ + grj^) + ... , 
y = vXcosel> + }vX^COS^<l>{i-i' + -ri'') 

+ ^vX'cos'>4>{5-iSf^ + t'+i4rf-5&v'^') •••• 

It will be observed that in the term in X* for x we have + ^rj'^ in the Cassini 

as against + 91?^ in the Gauss. 

The following value for the convergency on the Cassini is also found by 

continuing the method of Jordan : 

^ . , f X^ cos^ d> , 9V X* cos* <^ / ^2 2 2 v2\ 1 
■y = Xsin<^|l+ ^(l+rf)+ — -^(2-^^+4f-Sr)2/^) ...| . 

Now the work of J. I. Craig, p. 34, enables us to find the true bearing 
of the line on the spheroid to which a line, through any point on the 
projection, parallel to the axis of x corresponds. It is 

V cos 4) dyjd^) 
p dyjdX ' 
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Performing these differentiations we find 
dy\d<^ = - p\ sin ^ |i + ^^!f2!!i (2 - ^2 + 2^2) 

+ ^tE^t± (16 - 28/^ + t'^Z^i- - 6oy^e) ...} , 
120 ^ J 

dyldX^ + voos^L-^^^°t±t^-^^^^{^t^-t^+ir,^n..], 
I 2 24 J 

whence 

tan a = A sin ^ J I + ^ [i + r + rf)- 



X^ cos" 
15 



(2 +4^^ +2/" + 4^') ...I, 



. . , f k^ COS^ <b , », X" cos" <^ , ,, „ 2 ^2\ I 

a = A.sm<^ ji+ ^(i +v) + (2-i^ + AV - Sn i) ■■■] > 

agreeing with the convergency. 

When however we apply the same method to the Gauss Conformal 
Projection to find the true bearing on the spheroid to which a line of zero 
bearing on the projection corresponds we shall find 

. . , f \^ COS^ (^ , ,\ X" cos" 4> , _,;, o 2 J2\ 1 

a = X sm ^ j I + ^ (i + 377^) + — —^ {2-t^+ i5»?=- iS^;'^') ...\, 

differing from the value in the Cassini by the small quantity 

I rfX^ sin 4> cos^ ^ + 

This shews that in the Conformal System the lines of zero bearing do not 
correspond with the lines of zero bearing in the Cassini System, or with the 
lines of same bearing on the spheroid. The indicatrix, though a circle in 
the Conformal System, is at every point turned through the above small 
angle. This is a point not of practical importance perhaps, as the quantity 
in question is so small, but it has some theoretical interest in understanding 
the true nature of the convergency of map meridians. 

For the scale value along a line parallel to the principal meridian on 
the projection we also find from J. I. Craig's formulae 

fdx dy dx dy \ I f , fdyV ^ , , fdy^^ 
which in the case of the Cassini leads to 

%=0 = I + i'^' COS^ <i>{l+rf) + JiA." cos" <^ (5 - 4^2 + lOiJ^ - 2j^-fP) + ..., 

whereas in the Gauss Conformal we have 

^^=0 = ^ + ^-^'cos' <^ (i + yf) + ^X" cos" <^ (5 - 42'= + 141?^ _ 2?,yfP) + .... 
For any other value of fi we have approximately on the Cassini 
m =1 +^X^cos^</> (i +ij^)cos^y8+ .... 
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Some small errors in Helmert's Formulae 
for Convergency, etc. 

In checking the formulae for the convergency on the spheroid, etc. by 
Helmert's values the writer was led to discover some small errors in the 
latter, which should be pointed out, though they will hardly have any 
appreciable effect in practice. 

It will be found that the value for the convergency can be put into 
the form 

tan 7 = tan X sin <^ { I + ^r;^ cos^ <^\^ + tsV^ cos* </>X' (7 - 5?") + . . . } . 
Transferring this into Helmert's form, we write 

whence we have 

tan / = tan Z,.2 sin £,{i+ ^e^L^ cos^^^ [i + ^2 _ z^ (| - 1 cos^ B^)]}, 
as against 

tan /■= tan Zi_ 2 sin B^ {i + ^^ Z^ cos* B^ [i + ^e' (5-3 cos B^) 

-Z^(|-|cos^i?,)]}> 
as given by Helmert, p. 446, equation (7). 

The error seems to arise in the equations (2), (4), (5) and (6) on the 
previous page 445. In equation (2) the fourth line should read 

^Z*(3o-25 siii^i^cos^^- 1 6 cos* Z'), 
•instead of J^Z* (30 - 10 sin^ Fcos'F- cos^i^- 

Equations (4) and (5) are in error because the terms in ^ are omitted, 
though they are correctly given in equation (10), p. 442, for the relation 
between cos Z" and cosB^. These omissions affect equation (6) and the 
following equation which should read 

sin J82 = sin B^ {i-^e^ cos^ F[i + le" (1 + 3 sin^Z") + Z- sin^Z" 

+ i^Z^ sin^ Z'(i+ sin^ Z') - |Z* sin^ Z + Z* sin* Z" 
+ itf* (i + 2 sin^ F+ 5 sin* F)]}, 
tan / = - tan Z sin Zg {i + ^e^L^ cos* Z" [i + e^ + z^ (f - 1 cos^Z")]} . 
Helmert's investigation is extremely difificult to follow and the use of B^ 
the actual latitude of a point and F the latitude of the point where the 
perpendicular coordinate of the first point cuts the principal meridian, in 
the same formula, is very apt to lead to confusion. 

On the sphere the latitude of the foot of the perpendicular which 
determines the x coordinate in the Cassini Projection is given by 

tan F= tan Zj sec Z 
in Helmert's notation, or tan </>! = tan <^ sec A in ours. From this the well- 
known expansion for 1^1 - <^ is derived by trigonometry. For the spheroid 
the same formula can be put in the form 

tan ^i = tan <^ sec A{i +i7j^A^(i +^X^cps 2ij(>) ...}. 



7& EXPANSIONS 

In the same way the formula for the y coordinate can be put in the form 
sin {yfv^ = sin X cos <if> (i - xV'?^^* sin^ ^ ^^^^ '^)- 



List of Expansions used in these Papers. 

sin^ = ^-ie» + T-i^e=-^^e' + ^^^^e''-..., 
cos e= I - 1^^ + ^1^6^- ^e« + j^^^e»- ..., 
tane = e + |e^ + T\e= + ^Ve' + ^ff^e»+..., 

cote = |{i-ie^-J^e^-^^e«-...}, 

sec e = I + \&' + A^ + AV^' + ^VsV^' + - . 

sin-' e = e+ie> + l-^6'+^^6' + ^^1^6^+ ..., 

log, sin e = log e - 1^^ - xiTT^^ - WsT ^' - ftIw^' - • • • > 
log, tan e = log ^ + ^e^ + ^V^ + ^Ht ^' + xiwir^' + - ■ , 

\og,{i-e) = -{e+ie' + ^6'+i6* + ^6'+...}, 
ioge^g= 2 {e + i&'+^e' + ^^e^ + -...}, 

tan (6 + /^) = tan ^ + /% (i + t^) + /%V (i + /=) + J/%' (i + 4/^ + 3/*') + . . . , 
where i = tan ^. 
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